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ABSTRACT 


Search of an unknown space by a physical agent (such as an autonomous vehicle) is 
unique in search as the customarily most important goal (to reduce the computation time 
required to obtain the shortest distance) is not as important as minimal movement. There is 
a real-time aspect since the agent is actually moving; using energy each step of the way. 
Having limited energy resources and knowledge of the terrain (only adjacent nodes), the 
key factor for the physical agent’s search algorithm is reduction of steps. Hence, any 
heuristic that can help keep step count to a minimum must be considered. Korf and Shing 
addressed this issue in separate works. Both made use of known information about the 
frontier node’s distance from the current node in addition to a heuristic estimating the 
distance from goal. 

In this thesis, we present a simple genetics-based method to produce adaptive, efficient 
multi-heuristic search strategies for the real-time problem. Extensive empirical study 
shows that this approach produced search strategies with much better performance over 
existing search algorithms for most terrain types. The methodologies used to develope 
these improved strategies for our specific case, are also applicable to a multitude of real- 


time search/optimization problems in the general case. 
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I. INTRODUCTION 


Search of an unknown space by a physical agent (such as an autonomous vehicle) is 
unique in search as the customarily most important goal (to reduce the computation time 
required to obtain the shortest distance) 1s not as important as minimal movement. There is 
a real-time aspect since the agent is actually moving; having limited time to determine its 
next move and using energy each step of the way. This is in contrast to the traditional 
problem of search of known space for the shortest path which can be efficiently 
accomplished by A* search with a good heuristic estimating the distance to goal. The path 
is found without any movement. Although factors other than the actual distance from start 
and estimated distance from goal could reduce the number of nodes examined in the 
traditional problem, these factors usually increase the computational cost per node 
examined and produce paths that are longer than the shortest path which makes additional 
heuristics undesirable. Such is not the case for the real-time problem. 

The physical agent traversing a terrain in the real-time problem knows only its current 
position, the goal’s position, and whether adjacent and previously adjacent nodes are 
passable or not. It learns about the terrain only as it moves from node to node examining 
all nodes adjacent. Information about past nodes, visited or adjacent, can be stored to build 
up its knowledge base. Computational time to determine the next move is important, as 
stopping to compute before each move is undesirable. On the other hand, insufficient 
computations can result in unnecessary steps and wasted energy. 

Having limited energy resources and knowledge of the terrain the key factor in the 
physical agent’s search is the reduction of physical steps. In [Pa89], Papadimitrion and 
Yannakakis showed that the computational problem of deriving optimal search strategies 
for the real-tme problem is PSPACE-complete. Hence, any heuristic that can help keep 
step count to a minimum must be considered. Korf [Ko9Q0] studied this problem and 
developed the real-time-A* search, which uses the physical agent’s distance from the node 


(g (n) ) in addition to the distance from goal heuristic (4 (n) ) to determine the best next 


move by minimizing the objective function f(n) = g(n) +h(n) for every adjacent node 
n. Shing and Mayer [Sh91] developed persistence search which included a persistence 
factor (pf = 0 to 1) to bias the distance from current. The next move is determined by 
minimizing the objective function f(n) = pfxg(n) +h(n) for every frontier node n. 
Experimental results led to the conclusion that the pf could be adjusted to optimize search 
depending on terrain type and the density of obstacles. Details of these search strategies are 
in Chapter ITI. 

Extending on these works, we believe a combination of additional heuristics can be 
beneficial in minimizing physical agent steps. As the number of heuristics increases, it is 
essential to have some efficient means of assigning bias adjustments to various heuristics 
to optimize f(n) for different terrain types and densities of obstacles. If the combinatorial 
explosion required to produce all possible combinations of heuristics is not intractable, the 
required testing of each to select a best makes this means computationally prohibitive. 
Since enumeration is probably not possible, some random means of attaining the best 
combination seems to be the most plausible. DeJong [De75] made clear the advantages of 
genetic algorithms over purely random selection. 

In this thesis, we present a simple genetic algorithm based method to produce adaptive, 
efficient multi-heuristic search strategies for the real-time problem. Extensive empirical 
study showed that this approach produced search strategies with much better performance 
(reduced number of steps without prohibitive computation time) over existing search 
algorithms for most terrain types. The methodologies used to develope these improved 
strategies for our specific case, are also applicable to a multitude of real-time search/ 


optimization problems in the general case. 





Il. PROBLEM MODEL 


A. TERRAIN MATRIX 

To best demonstrate the effectiveness of the multi-heuristic search strategies produced 
by a genetic algorithm, we chose to apply the strategies to random obstacle distributions in 
the form of a two-dimensional 64x64 grid of squares (nodes). Nodes can be either free or 
Piriacles: movement can be in eight directions through free spaces only. A perimeter 


surrounding this grid is a solid row/column of obstacles. The distance from a node to its 


horizontal/vertical neighbor is 1.0; to its diagonal neighbor is J/2. The total distance 
traveled from start to goal according to any search scheme is the sum of each of these 
individual steps. The effectiveness (fitness) of a specific search scheme is the ratio of the 
shortest path length from start to goal divided by the distance traveled. Given as a 
percentage, 100 is the best possible; meaning the distance traveled is equivalent to the 
shortest path. Specific nodes of the grid are be identified by Cartesian coordinates with the 
left border column being the y axis and the bottom border row being the x axis. The lowest 


left node 1s (1,1); the top nght is (64,64). 


B. DENSITY MATRIX 

The 64 by 64 search space grid is divided into 16x16 density blocks, each containing 
4x4 nodes and having a specified block density. Block densities range from 0-15. A block 
density of 9 means that, on average, nine of the block’s 16 nodes will be an obstacle (chosen 
at random). These density blocks are numbered from (0,0) to (15,15) where (0,0) is the 
bottom left and (15,15) is the top right. Start and Goal positions are specified by density 
blocks. Most of the block density distributions used will have a start block of (2,2) and a 
goal block of (13,13). The specific start/goal node is located randomly in that block. See 
Appendix A. 


C. TERRAINS USED 


There are ten different density distributions that were used for training and testing. The 
block densities, once set, remain unchanged from the start of training through testing. 
Although the block densities remain constant, actual obstacle placement is determined 
stochastically.and changes from-run to run. The point is to investigate the adaptability of 
genetic algorithm to produce the best strategy to direct the search through terrains where 
the general density distribution is known but actual obstacle placement is not. The first six 
terrains are considered natural terrains since they closely resemble actual topological 
conditions. The start density block is always (2,2) unless otherwise stated. The goal density 


block is always (13,13) unless otherwise stated. 


1. Central mountain 
The highest density, 15 (denoted as f in Figure 1), is in the center with a gradual 
decrease towards the lowest density, 1, on the outer edge. Figure 1 shows the density 
distribution of the terrain in hexadecimal. Transit from start to goal requires a search 


scheme to find the most efficient way around the mountain. 
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Figure 1: Central Mountain Terrain 








2. Single Left Ridge 
This terrain has a high density (15) ridge starting from left center moving 
horizontally out past the grid’s midpoint. There is a gradual decrease in density down to 2 
as the distance increases from ridge center (Figure 2). Transit through the ridge is not 


possible. 
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Figure 2: Single Left Ridge Terrain 


3. Single Right Ridge 
This terrain is similar to the Single Left Ridge but in the opposite direction 
(Figure 3). This is a much more difficult situation since the physical agent must move away 


from the goal to find its best route. 


NN GFW SOWW O BOUT AN 


NN FwWoO SOWH OO GOW ITAA 
NN Pwo SOW H OO BOM GAA 
NN FWO BTOWH O BOW ETNA 
NN PFW BPOWH OO BDTOU TNA 
NN Fw POW O BOM TAN 
NN FW BGOWH O BOW ENA 
NN PFWO BGOWH O BOWIE NAN 
NN FWO TOW O TOM TAN 
NN PFWO PTOWH OO TOW TAN 
NNPFWO BSOWH O BOM TAN 
NNN TWO BOO BMW YENNN 
NAUNNFWO BBDWENNN 
NNNNN FWDWOWIENNANN GO 
NAUNNNN FWOWTNNNANNN GS 
NAUNNANNN ES ENN 


Figure 3: Single Right Ridge 


Double Ridge 


4. 


This terrain has density areas producing a right ridge on top of a left ridge. Ridge 


densities are 15 with a valley of 2 in between (Figure 4). An s-shaped path to get from start 


to goal is required to transit this terrain. 
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Figure 4 : Double Ridge Terrain 





5. Single Left Plateau 
This terrain is characterized by a large area of high density (10) (dense but 
passable) starting from left center moving horizontally out past the grids midpoint (Figure 
5). The start/goal density blocks are (4,0)/(11,15). A successful transit can consist of either 


direct passage through the plateau or circumnavigate. 
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Figure 5: Single Left Plateau Terrain 


6. Single Left Plateau with Ridges 
The same as single left plateau, except it has a ridge of higher density (12) (hard 
to pass) along the plateaus perimeter (Figure 6). The start/goal density areas are (4,0)/ 


(11,15). Circumnavigation will usually be the only viable option. 
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Figure 6: Single Left Plateau with Ridges Terrain 


7. Random Terrains 
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from 0 to 15 with equal probability. Shown in Appendix B (Figures 30, 31, 32 & 33), these 


grids were used to test the effectiveness of the different search strategies through unnatural 


terrains. 





Il. BACKGROUND WITH DISCUSSION PERTINENT TO MODEL 


A. KNOWN TERRAIN SEARCH 


A guarantee of optimal path from start to goal is usually the main concern in known 


terrain search. 


1. A* Search 
[Ha68] - Widely accepted as the best algorithm for finding the shortest path in a 
known search space, it uses the actual distance from start and a heuristic estimating the 
distance from goal. The object of any frontier node is given by the following equation: 
f(n) = s(n) +h(n) (Eq 3.1) 
where s(n) = the actual distance from start to n, the frontier node, and h(n) = the 


heuristic estimated distance from n to goal. Although guaranteed to find the optimal path, 
assuming that the heuristic estimate is always less than the actual, it 1s not required to 
minimize the number of nodes examined. Using Euclidean distance as the distance to goal 
heuristic, A* search provided the shortest path for each terrain used in our experiment to 


compare the effectiveness of each search scheme. 


B. UNKNOWN TERRAIN SEARCH BY A PHYSICAL AGENT 

Assuming limited sensory range, the physical agent cannot find the shortest path 
without excessive moves. Although the shortest path would be nice, more important to the 
search schemes success is the energy expended/time spent finding a satisfactory path. 
Distance traveled is of major concern as is time to determine next move (related to this is 
the computational limitations of the physical agent). The following algorithms have been 
considered in an attempt to find the best. In all equations, n represents one of the frontier 


nodes on the frontier list unless otherwise stated. 


1. Hill-climb Search 
[Wi92] - Depth first search with each move determined by the best (least distance 
from goal) adjacent frontier node; n In the following equation: 
f(n) = h(n) (Eq 3.2) 
When no frontier nodes are adjacent to the current, the algorithm backtracks until adjacent 
frontier nodes are encountered. This search scheme minimizes jumps in search of the best 
frontier node, but pays the price in extra steps due to unchecked persistence on initially 


good but eventually poor tracks and the often resultant steps required to backtrack. 


2. Real-time-A* Search 


[Ko90] - Uses distance from current (actual) and distance from goal (heuristic) to 


determine best next move. This search only looks at adjacent nodes (frontier and visited). 
In the following equation n stands for adjacent non obstacle nodes. 

fn) = a(n) +h(n) (Eq 3.3) 
g (n) 1s the actual distance from the current to the adjacent node n. The h(n) is a heuristic 
predicting the distance from n to the goal. Initially, h (7) is calculated by using Euclidean 
distance in our example. The algorithm picks the adjacent node with the best f(n) . Before 
moving, the value of the f(n) of the second best adjacent node is stored in the current node. 
This stored value will, in future computations, be the node’s h(n). This value remains 
constant until the node is visited again. This well conceived search scheme requires 
minimal computations and memory, but is subject to wasted moves when drawn into local 


traps. 


3. Best-first Search 


[Wi92] (modified for physical agent) - Uses only the distance from goal heuristic 


(Euclidean distance) to select the next move. 


f(n) = h(n) (Eq 3.4) 
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Once selected, it uses the shortest path through visited nodes (known search space) to 
travel the distance from the current node to the selected frontier node. Although, after each 
move it is at the best known location, the cost of getting there can be expensive. In worst 
Case situations it can end up jumping large distances back and forth while zeroing in on the 


goal. 


4. Persistence Search 

[Sh91] - Similar in concept to Real-time-A*, it uses the distance from goal 
heuristic and a weighted distance from current to determine its next move. Unlike Real- 
time-A*, it makes more use of known information; resulting in better moves, but decreased 
computational efficiency. The distance from goal is Euclidean. The distance from current 
to frontier is the shortest path through visited nodes as in Best-first, but this distance is 
weighted and used in determining the next move. The object of a frontier node 7 is given by 
f(n) = pfxg(n) +h(n) (Eq 3.5) 
where g (n) =shortest distance from current position to n through visited nodes, h(n) = 
Euclidean distance from n to goal. A persistence factor (pf = 0.0 to 1.0) is added to vary 

the relative contribution of each of the heuristics to the determination of next move. 
Distance from current, assumed to be always less pertinent, can be reduced in importance 


in comparison to distance from goal. 


C. GENETIC ALGORITHMS 

Genetic algorithms, developed by John Holland [Go89] and his associates, are based 
on the laws of natural selection and survival of the fittest. Subjecting a population (animals, 
search schemes, etc.) to environments where fitness for survival is required, individuals 
best suited for survival will flourish and reproduce while individuals lacking the diversity 
required to continue in all possible environments will discontinue. 

The key to the success of a population is its robustness [Go89]. An individual, and 
therefore a population, is made up of traits which are derived from specific genes in the 


individuals chromosome [St77]. Applicable traits in the animal world are weight, height, 
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leg length, neck length, etc. A combination of these traits describe an individual. Extremes 
in any one trait usually means more specification and added survivability in a limited range 
of environments, whereas moderation in traits means added adaptability for diverse 
environments. The key is to find the balance of these two in a population to give it proper 
robustness. Example: the giraffe can afford to be specifically designed for reaching (long 
neck and legs) because it doesn’t face diversities in environment that would require escape 
through low canopy jungles. It is perfectly adapted for life on the plains with occasional 
trees. 

Similarly, search strategies can be very specialized in simple environments. Search 
through a low density (of obstacles) terrain can be successfully accomplished with 
efficiency and consistency by a simple hill climb algorithm (only one trait, distance to goal 
of adjacent nodes, is important). Search problems involving more complicated and 
diversified solutions require the proper balance of traits (heuristics) to solve. Simple direct 
“hill climbing” approaches can often result in searching locally optimal blind alleys. 

One possible means of developing the balance of traits required to avoid getting stuck 
in the local minimum is to enumerate all possible combinations. This would most assuredly 
find the optimum, but in many problems the combinatorial explosion of possibilities make 
this method prohibitive. Purely random combination of trials is a possibility that seems to 
avoid both the local minima and the combinatorial explosion problems. But on further 
examination, it suffers the same drawbacks as enumeration, in that there are only a limited 
number of trials possible whether you look at them in order or at random. Genetic 
algorithms use randomness as a tool in a direct search for the optima. Promising potential 
solutions can be searched in parallel while feedback information 1s used to select the next | 
partially random strategy. The results, as evaluated by DeJong [De75], show the superiority | 
of genetic algorithms over purely random. } 

The basic genetic algorithm makes use of a population of individuals (usually binary | 
strings of fixed length) that are made up of the traits pertinent to the problem (traits are 


usually represented by a fixed number of the bits in specific locations). Three genetic 
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operators are used to transform the original (randomly generated) population into an 
optimal one: selection, crossover, and mutation. 

The Fitness of an individual of the population is established by some form of 
evaluation function. One scheme is to compare each to a known optimum, assigning higher 
fitness to ones approaching the optima. This evaluation can also be averaged over some set 
number of trials (cycles) for each individual and then assigned as the fitness before forming 
the next generation. 

Each new generation of the population is formed by stochastically selecting 
individuals from the prior population. Higher fitness individuals have a higher chance of 
being selected. Reproduction is performed by randomly pairing selected individuals for 
crossover and mutation. 

Crossover is performed at a random point in the binary string. The two selected strings 
interchange their tail sections at the crossover point to from two new individuals. The 
crossover point can be anywhere from 0 to the last bit. For example, let the two selected 
strings be QOOOO000 and 11111111, and let 5 be the crossover point. Then the crossover 
operation will produce the new strings 00000111 and 11111000. In general, crossover 
forms two new individuals with one hopefully having all the best from its two parents. 

Mutation is a bit by bit operator that takes each individual and randomly (with a 
specified probability) decides if each bit will be changed? For example, changing the 
second bit of the string 00000111 by mutation will result in the new string 01000111. This 
genetic operator, as in nature, ensures that populations maintain adaptability even when 
specialization is the rule. An extremely high mutation probability regresses the genetic 
algorithms to a uniform randomly distributed population, a very low one reduces the 
populations adaptability. A happy medium seems to be in the 0.01 to 0.001 range for 
probability of individual bit mutation. 

A myriad of variations are possible to improve the performance and robustness of the 
genetic algorithms. For the purpose of our research, these were found to be unnecessary, 


and will not be covered in this discussion. 
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IV. FACTORS RELEVANT TO SEARCH 


A. DISTANCE FROM START 

This is usually the actual shortest path from the start node to the considered frontier. 
Currently believed to be useless in a real-time environment, it should be selectively 
eliminated by natural selection as the genetic algorithm trains. For our implementation, it 
is approximated by computing the Euclidean distance from start to frontier. It may be 


significant in some of the more complex terrains that require a switch back. 


B. DISTANCE FROM CURRENT 


The distance from the current node to the frontier node; important in Real-Time-A* 
and Persistence Search to determine if backtracking is worth the steps required. It is the 
actual distance computed as the actual steps required to move from the current node to the 


frontier. 


C. DISTANCE FROM GOAL 


The Euclidean distance from the current node to the goal node. This heuristic is usually 
considered important in any search. It is used in combination with “distance from current” 


for Persistence Search, and by itself for Best-first Search. 


D. CROWDING 


The crowding parameters, crowding sides and crowding diagonals, are an attempt to 
assist the physical object in avoiding areas of increased obstacle density. This reduces 
exploration of paths through high density areas, favoring the safer path of increased options 
available in the open space. The parameters are separated in case one is more appropriate 
than the other. Both would be much more effective without the self imposed constraint of 
physical object perception only being adjacent nodes. If all nodes adjacent to the frontier 


node could be seen, these factors importance would increase significantly. 
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1. Crowding sides 
This heuristic examines the frontier node’s known horizontal/vertical neighbors 
to count the number of obstacles. Nodes with more known obstacle neighbors are less 


desirable. The minimum value is 0 and 4 is the maximum. 


2. Crowding diagonals 


This 1s similar to the previous parameter with the count being made of the frontier 


node’s diagonal neighbors. 


E. MOVE AWAY FACTOR 


It attempts to continually reduce the search space by reducing desirability of nodes 
that increase the x and/or y difference between the current and goal nodes. Increasing the x 
or y distance counts as 2, increasing both counts as 4, and no increase results 1n the heuristic 


having a value of 0. 


F. MOMENTUM 

This heuristic attempts to avoid zigzag by making forward (in relation to last move) 
nodes the most desirable. It should be useful in valley/ridge terrains where the best path is 
straight through the valley. By maintaining momentum, the physical object avoids steps 
wasted in popping in and out of each crevice which has nodes closer to the goal. Straight 
ahead movement results in a value of 0, a 45° shift makes it 1, a 90° shift is 2, and a 135° 
shift or non-adjacent move results in a value of 3 (making only the adjacent nodes subject 


to change after a move). 
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V. PROGRAM DEVELOPMENT 


A. DATA STRUCTURES 


1. Node structure 
The 64x64 grid is internally represented as a 66x66 two dimensional array (the 
perimeter nodes are all marked as obstacles) made up of pointers to node records. The 
records store information pertinent to terrain, search (heuristics), graphic display, and 
pointers to other node records (used in the program for various dynamic structures). The 
heuristic values stored include distance from start, distance from goal, distance from 
current, crowding sides, crowding diagonals and subtotal. No other node records are used 


in the program; other structures requiring nodes are set up using pointers to these records. 


2. Population structure 
A 32 member array of individual records makes up the population. Each stores 
the individual’s fitness and its chromosome which contains biases for each search 
parameter. The chromosome is a 32 bit unsigned integer; subdivided into eight four-bit 


unsigned integers, it holds up to eight heuristic bias factors with a range from 0 to 15. 


XXXXXKXX 


place holder small 
start distance factor 


goal distance factor 
current distance factor 


crowd sides factor 


crowd diagonals factor 
move away factor 


momentum factor 
Figure 7: Individual Chromosome Structure 
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3. Frontier Heap 
Implemented as an array of pointers to node records. Functions to manipulate this 


min-heap are in theap.c. 


4. Mate Heap 
An array of 31 integers. The first is a count of the heaps members, the other 30 
are chosen at random in the range of 0 to the sum of all individual fitnesses. They are used 
to randomly select individuals for reproduction. Functions to manipulate this heap are in 


eheap.c. 


B. SEARCH ALGORITHMS 

Algorithms used for this analysis, with appropriate modifications, are covered in this 
section. Each node, represented as a record, has a number of fields that are used to store 
needed information. When a function operating on a specific node v is used, a read or write 
to the appropniate field takes place. For example: in the A* algorithm the following 
calculation takes place; f(v) = s(v) + h(v). The value s(v) calculated earlier was stored in 
v.s (the nodes dist_from_start field), h(v) is calculated and stored in v.h (the nodes 
dist_from_goal field), and f(v) is stored in v.f (the nodes subtotal field). 


1. A* Search 
This search (Figure 8) will find the shortest path from start to goal if a path exists. 
The heuristic estimating distance to goal (h(v)) is a lower bound of the actual cost of the © 


optimal path from v to the goal. 


2. Hill-Climb Search 
This search (Figure 9) always moves forward until there is no where to go. It then 
back-tracks the way it came until a move is possible. It is a depth-first search with a 


heuristic to determine the best move to advance it to the goal. 
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a_star_search 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
ay) 
(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 


current := start 
while current != goal do 
for all nodes, v, adjacent to current do 
if UNTOUCHED 
s(v) = current.s + distance(current, v) /* Euclidean */ 
f(v) = s(v) + h(v) /* h(v) is the Euclidean dist to goal */ 
add to frontier heap 
elsif FRONTIER 
if s(v) > current.s + distance(current, v) 
update s(v) and f(v) 
update position in frontier heap 
endif 
endif 
if frontier_heap is empty 
return BIG_NUMBER /* there is no path from start to goal */ 
endif 
end for loop 
current := top(frontier_heap) 
end while loop 


retum goals 


Figure 8: A* Search Algorithm 
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hill_ climb_search 
(1) current := start 


(2) while current != goal do 


(3) best := dummy_node /* f(dummy_node) = BIG_NUMBER */ 
(4) for all nodes, v, adjacent to current do 

(5) if UNTOUCHED 

(6) f(v) := distance(v, goal) /* the Euclidean dist to goal */ 
(7) mark v as FRONTIER 

(8) endif 

(9) if FRONTIER and f(best) > f(v) 

(10) best :=v 

(11) endif 

(12) end for loop 

(13) previous_current := current 

(14) if best != dummy_node 

(15) current := best 

(16) curlrent.predicessor := previous_current 

(17) elsif current != start 

(18) current := current.predicessor 

(19) else 

(20) retum BIG_NUMBER 

(21) endif 

(22) current.dist_traveled := 


previous_current.dist_traveled + distance(current, previous_current) 
(23) end while loop 
(24) return goal.dist_traveled 


Figure 9 : Hill Climb Search Algonthm 
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3. Real-Time-A* Search 
This algorithm (Figure 10) is in accordance with Korf’s description [Ko90]. For 
our implementation, the node array was used to store the h value since it was already in 


place, negating the necessity for a hash table. 


4. Best-First Search 
This search (Figure 11) always goes to the best (minimum h(v)) node regardless 
of its distance from the current node. It is possible to implement as a specific case of the 


multi-heuristic search (Figure 13). 


5. Persistence Search 

Shown in Figure 12, gf + hf are intended to effectively replace/descretize/expand 
the persistence factor, pf, in the original work ([Sh91] equation 3.5). pf can have any value 
between 0.0 and 1.0. We found that an infinite range of possibilities for this factor was not 
required. A descrete, yet sufficient, span can be obtained by setting gf and Af to any 
number of possibilities where gf < hf. Setting hf to 15 and incrementing gf from 0 to 15 

gives us the equivalent of a 0.0 to 1.0 range incriminated by 0.067. 
f(v) = efxg(v) +hfxh(y) (Eq 5.1) 
There is also now the expanded capability of having the g(v) be the more 
important factor in the search (gf> hf). This search can also be implemented as a specific 


case of the multi-heuristic search. 
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real time_astar_search 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 


(8) 

(9) 

(10) 
(11) 
(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 


(22) 
(23) 


current := start 
best := dummy_node 
second_best := dummy_node /* f{(dummy_node) = BIG_LNUMBER */ 
while current != goal do 
for all nodes, v, adjacent to current do 
if UNTOUCHED 
h(v) := distance(v, goal) /* Euclidean */ 
/* else h(v) is already set */ 
endif 
g(v) := distance(current, v) 
f(v) := g(v) + h(v) 
if best.f > f(v) 
second_best := best 
best :=v 
elsif second_best.f > f(v) 
second_best :=v 
endif 
end for loop 
previous_current := current 
current := best 
previous_current.h := second_best.f 
current.dist_traveled := 
previous_current.dist_traveled + distance(current, previous_current) 
end while loop 


retum goal.dist_traveled 


Figure 10: Real Time A* Search Algorithm 
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best_first_search 
(1) current := start 


(2) while current != goal do 


(3) for all nodes adjacent to current do 

(4) if UNTOUCHED 

(5) h(v) := distance(node, goal) /* Euclidean */ 
(6) add to frontier_heap 

(7) endif 

(8) end for loop 

(9) if empty(frontier_heap) 

(10) retun BIG_NUMBER /* no solution */ 
(11) endif 

(12) v.dist_traveled := current.dist_traveled + g(v) 


where g(v) is the shortest distance through known paths from 
current to frontier node. 

(13) previous_current := current 

(14) current := top(frontier_heap) /* minimum h(v) */ 

(15) end while loop 

(16) return goal.dist_traveled 


Figure 11: Best First Search 
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persistence_search 


(1) 
(2) 
(3) 
(4) 


(5). 


(6) 
(7) 
(8) 
(9) 
(10) 
(11) 
(12) 


(13) 
(14) 
(15) 
(16) 
(17) 


Current := start 
while current != goal do 
for all nodes adjacent to current do 
if UNTOUCHED 
h(v) := distance(node, goal) /* Euclidean */ 
add to frontier_heap 
endif 
end for loop 
if empty(frontier_heap) 
return BIG_NUMBER  /* no solution */ 
endif 
find the frontier node, v, that minimizes the equation: 
f(v) := gf * g(v) + hf * h(v) where g(v) is the shortest distance through 
known paths from current to frontier node. gf and hf, set before search, 
are bias factors used to vary the relative importance of g(v) and h(v). 
They can have a value from 0 to 15. 
v.dist_traveled := current.dist_traveled + g(v) 
current:= v 
remove current from frontier_heap and update 
end while loop 


return goal.dist_traveled 


Figure 12: Persistence Search Algorithm 
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6. Multi_heuristic Search 
This is the general algorithm (Figure 13) enstantiated in our case to handle five 
stable_heuristics and two unstable_heuristics. Stable_heuristics being ones that have 
values that will not change if more than two steps away from the current node. They include 
Euclidean distance from goal (hg), Euclidean distance from start (hs), crowd sides (hes), 
crowd diagonals (hcd), and momentum (hm). The subtotal fs (v) is calculated using these 
functions multiplied by their respective bias factor and stored in v.subtotal. 
fs(v) = hefxhg(v) +hsfxhs(v) + hesfxhes(v) +hedfxhcd(v) +hmfxhm(v) (Eq 5.2) 
Unstable_heuristics have values that are liable to change as the current node 
changes. Examples in our case: distance from current (hdc) and move away (hma). The 
algorithm minimizes equation 5.3 using the efficient “branch-and-bound” search through 
known (visited) nodes described in section 4.3 of [Sh91]. 
f(v) = fs (v) + hdcf x hdc (v) + hmaf x hma (v) (Eq 5.3) 
The hsf, hef, hdcf, hcsf, hcdf, hmaf, and hmf are bias factors that correspond with the 
individual chromosome’s lower 28 bits which are set during training. The highest four bits 
are, in our implementation, a place holder for future additional heuristics since only seven 
applicable heuristics were identified. Note that the Best-first and Persistence Search could 
be implemented as special cases of the multi-heuristic search algorithm. Best-first uses an 
individual chromosome input of 00100000 (the third factor being Agf). Persistence Search 
uses an individual chromosome input of OOxyO000 with x and y varying from 0 to 15 
(fourth factor being hdcf). 
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multi_heuristic_search 


(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 
(11) 
(12) 
(13) 
(14) 
(15) 
(16) 
oo) 


(18) 


(19) 


(20) 
(21) 


Current := start 
while current != goal do 
for all nodes v within 2 moves of current do 
if adjacent and UNTOUCHED 
v.subtotal := inner_product(stable_heuristics * respective_biases) 
add v to frontier_heap /* min subtotal node on top */ 
elsif FRONTIER 
if any stable_heunistics of v have changed 
v.subtotal := v.subtotal + adjustment 
update position in frontier_heap 
endif 
end if 
end for loop 
if empty (frontier_heap) 
return BIG NUMBER §_/* no solution */ 
endif 
find frontier node, v, that minimizes 
f(v) = v.subtotal + inner_product(unstable_heuristics * respective_biases) 
v.dist_traveled := current.dist_traveled + g(v) 
where g(v) is the shortest distance through known paths from current to 
frontier node. 
current := v /* and remove v from heap */ 
end while loop 


return goal.dist_traveled 


Figure 13 : Multi Heuristic Search Algorithm 
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C. GENETIC ALGORITHM 


The task of the genetic algorithm is to find the combination of the seven bias factors — 
that will result in the optimum search scheme. The values of these seven bias factors are 
stored in a single individuals chromosome. Application of genetic operators to a population 
(32 in our case) of these individuals will, after numerous iterations, produce our desired 
optimal individual. 

The genetic algorithm, described in this section, is invoked during training after some 
predetermined number of cycles (making up one generation). The input population will 
have a fitness value (ability to get through the terrain) assigned to each of it’s 32 individuals 
(details of this process are described in the next chapter). This fitness value and the 
individual’s chromosomal make-up are required by the genetic algorithm. 

Our algorithm (Figure 16) makes use of the three genetic operators: selection, 
crossover, and mutation. The implementation is similar to the algorithm presented in 
chapter one of the text by Goldberg, [Go89], with the additions of allowing the best two 
individuals to go unchanged and an average of one out of seven of the remaining not going 
through crossover. The result is similar to De Jong’s R3 elitist model [De75]. Examples of 


Our crossover implementation are detailed in Figures 14 & 15. 


Alleles are represented in hexidecimal 


Before allele crossover: 55555555 / 88888888 
Randomly picked crossover allele position is 3 (4th allele) 
After allele crossover: 55558888 / 88885555 


Figure 14: Allele Crossover Example 


The 4th allele is expanded out into binary representation 


Before bit crossover: 555 0101 8888 / 888 1000 5555 
Randomly picked crossover position between bits is 2 
After bit crossover: 555 1001 8888 / 888 0100 5555 


Figure 15 : Bit Crossover Example 
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genetic_algorithm 


(input is a population of individuals) 


(1) 
(2) 
(3) 
(4) 


(5) 


(6) 
(7) 


(8) 

(9) 

(10) 
(11) 
(12) 
(13) 
(14) 
(15) 
(16) 
(17) 
(18) 
(19) 
(20) 
(21) 
(22) 


total_fitness := all individual fitnesses added together 
select 32 individuals as follows /* selection */ 
best := individual with the highest fitness 
second_best := individual with the second highest fitness 
/*second_best must be distinct from best */ 
stochastically select 30 individuals with higher fitness individuals having 
the greatest chance of selection 
end selection 
create new_population with these 32 individuals 
pair individuals in such a way that it is unlikely that an individual is paired 
with itself; pair best with second_best 
for each individual pair, except best and second best, do 
randomly pick crossover allele position /* crossover */ 
ifnotQ /* 0 means no crossover */ | 
exchange all alleles after the crossover allele 
randomly pick crossover position between bits of selected allele 
ifnotOor4 /* 0 or 4 means crossover does not breakup the allele */ 
exchange bits after the crossover position between bits 
endif 
endif 
for each gene of the individual /* mutation */ 
invert bit if random < prob of mutate 
end for loop 
end for loop 
add best + second_best to new_population as individuals 0 & 1 respectively. 


retum new_population 


Figure 16: Genetic Algonthm 
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VI. TRAINING 


Training of the population is analogous to selectively breeding a random group of 
asexual organisms to obtain superior capability in a specific area. The capability you wish 
to optimize is transit from start to goal in the least number of steps. The specific area is a 
specific terrain layout where you have an idea about general areas of density, but have no 
information about the location of specific obstacles. 

The first step is to generate a series of specific terrains from your general idea of the 
densities. This can be done by placing obstacles in each area if a randomly generated 
number is less than the specified density. In our implementation, we simply loop through 
the 64x64 node array assigning each nodes state to OBSTACLE if the random number is 
less than the density value of the corresponding density block. The second step is to 
generate a population of 32 individuals giving them randomly generated chromosomes. 
Now the training begins (Figure 17). In all our work, we used 1000 generations with five 
cycles (trials) per generation. 

The returned population evolves through the trials of 5000 terrains. One of the 
individuals of this population is likely to have a chromosome that approximates the 
optimum combination of bias factors. Identification of this individual is accomplished 


during testing. 
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training 


(1) for the number of generations do 


(2) for the number of cycles do 

(3) loop until a successful A* search 

(4) create a terrain from the density_array 

(5) shortest_path := A* search 

(6) end until loop 

(7) run each individual through the terrain accumulating its fitness_sum by 


comparing its actual path to the shortest path 


(8) end for loop 

(9) compute each individual’s average fitness from fitness_sum and 
number of cycles 

(10) apply the genetic algorithm to the population 


(11) — end for loop 


(12) return a trained population. 


Figure 17: Training Algorithm 
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VII. TESTING 


Testing of the trained populations was performed by comparing the search conducted 
by the best individual in each population to searches accomplished using Hill-climbing, 
Best-first, Real-time-A*, and Persistence search. The following equation was used to 
compute fitness for all search schemes: 

fitness = integer (((Shortestpath) + (actualpath)) x 100) (Eq 7.1) 
Each search scheme was tested on 500 distinct terrains produced using the corresponding 
density matrix. 

Before testing, the best of each population was chosen by running the population 
through 50 distinct terrains. The individual with the highest fitness was chosen to represent 
the GA-trained population. The best values for distance from goal and distance from 
current bias factors for the Persistence search were determined by running 32 combinations 
(chromosomes of 0000000 to OOffO000 and OOO0f0000 to OOffOOOO) through 50 distinct 
terrains. From this, the best combinations of the two factors was used to represent 
Persistence search. 

The GA-produced best individual, Persistence best, Hill-climb, Real-time A*, and 
Best-first schemes were then all used to find a path in the 500 separate terrains. Average 
fitnesses over the 500 were assigned and a comparison of these fitnesses is presented in the 


results. 
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VIll. EXPERIMENTAL RESULTS 


The fitness of each search scheme in these results is the number of its required steps 
divided by the minimum steps possible, averaged over the 500 terrains used for testing. 
Fitness is presented as a percentage, with a 100% search scheme being one that can, on the 
average, search a terrain type in the minimum steps possible. In general, the easier the 


density layout of the terrain, the higher the fitness will be. 


A. NATURAL TERRAINS 


A graph comparing the fitness of applicable search schemes is presented for each 
natural terrain density layout (Figures 18 - 23). The following discussion is pertinent to 


each of these comparisons. 


1. Central Mountain 
This graph (Figure 18) shows that this terrain is only moderately hard for all the 
search schemes. Persistence search with a distance to current factor (gf) of 15 and a distance 
to goal factor (hf) of 11 (gf/hf = 15/11) was the best of the conventional search methods. 
The genetic algorithm produced an individual with chromosomal make-up of f00732b9 
(see figure 7, page 16 for breakdown) which performed 1.20 times better than the best 
conventional. Driven more to the goal by the move-away heuristic than distance to goal, 


this scheme was better equipped to avoid the congestion of the central mass. 


2. Single Left Ridge 
Overall this terrain was a little harder than the Central Mountain but was still 
handled moderately well by all search strategies (Figure 19). The best conventional was 
again persistence search using a gf/hf ratio of 15/6. The genetic algorithm scheme 
(f00828ff) had a fitness 1.16 times as good as the best Persistence and 1.28 times better than 


the next competitor (Hill-climbing). 
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Best First 
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Figure 18 : Central Mountain Results 
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Figure 19: Single Left Ridge Results 
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3. Single Right Ridge 

As could be predicted, this was a hard problem (Figure 20). The normal search 
schemes tend to spend a lot of time searching every possible route that was most direct to 
the goal. They would get stuck under the ridge with no way around, except back the same 
way they came. The best of these was Hill-climbing since it probably doesn’t waste a lot of 
steps backtracking. The best Persistence search had a gf/hf ratio of 2/15 showing it’s 
favoritism for a no backtrack approach. The search scheme produced by the genetic 
algorithm was superior to all by a multiplication factor of 1.26. Its chromosomal make-up 
was f00c2ca8. This scheme considers distance to goal to be not significant. It instead uses 
move away factor as the drive toward the goal. As the Persistence search with its gf/hf ratio 
of 2/15 the genetic algorithm produced scheme considered the amount of backtracking a 


major factor. 


Best First 
Persistence 
Hill Climb 
Real-time-A* | 


GA-produced | “a 
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Figure 20: Single Right Ridge Results 
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4. Double Ridge 

This terrain was a very difficult problem (Figure 21). Requiring negotiation 
around two ridges which involved a switch-back away from the goal, none of the search 
schemes were over 20% fitness. The average of the five schemes was 10%. The best 
Persistence, with a gh/hf ratio of 6/15, was roughly equivalent to Hill-climbing. The genetic 
algorithm generated scheme with a chromosomal make-up of f83b19bc was the best 
strategy with a fitness 1.21 times better than Persistence. Here is an example where distance 
to start was of significance; probably helpful in influencing the search to make the switch- 
back away from the goal. Move away factor was a major influence In striving toward the 
goal, backtracking was determined to be non-productive, but maintaining momentum was 
found to be important. It’s interesting to note that diagonal crowding was considered more 
important than side crowding (no explanation). The complexity of this scheme with the 
subtle interaction between these differing bias factors helps to confirm the necessity of a 


genetic algorithm to sort them out. 


Best First 
Persistence 
Hill Climb 
Real-time-A* 


GA-produced 





20 30 40 50 60 £470 80 90 100 
FITNESS (percent) 


Figure 21 : Double Right Left Ridge Results 
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5. Single Left Plateau 

This terrain was slightly easier for all the search strategies (Figure 22) although it 
presented a unique problem. The through the plateau route is possible but requires 
numerous explorations. Circumnavigating the plateau saves exploration steps but costs in 
the distance required. Since each of the 500 terrains had varying obstacle placement, we 
suspect sometimes it was best to transit through and other times better to go around. Since 
no general path was consistently optimal, the genetic algorithm had to develope a scheme 
that was equally effective for both routes or concentrate on perfecting one. In either case, 
its performance was again superior by a significant margin (multiplicative factor of 1.19). 
The resultant chromosomal make-up was f05e884f. The next best was Persistence with a 


gf/hf ratio of 11/15. 
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Figure 22: Single Left Plateau Results 
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6. Single Left Plateau With Ridges 
This terrain adds topological characteristics that favor circumnavigation as a 
search strategy. The genetic algorithm produced scheme, with a chromosomal make-up of 
f07c033d, was the best by a multiplication factor of 1.17 over the next best competitor 
(Figure 23). Momentum being the most important factor, it probably helped keep the search 
moving horizontally until clear of the plateau. Persistence was again the second best with 


a gh/hf rato of 11/15. 
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Figure 23 : Single Left Plateau With Ridges Results 


7. General Comment 
The genetic algorithm was extremely successful in producing the best search 


strategies for all natural terrains. 
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B. RANDOM TERRAINS 


Although the genetic algorithm produced search schemes where clearly superior for 


the natural terrains, we wanted to test their viability on randomly generated terrains. 


1. Random One / Random Two / Random Three 


The results from these three terrains showed that the search heuristics produced 
by genetic algorithms was of minimal value (Figures 24 to 24). These were all simple 
problems with the average fitness for all the search schemes being 64%. Fitness varied little 
between search strategies with a maximum of a 8% difference between the best and the 
worst. The genetic algorithm produced scheme was 1.02 (Random One), 1.01 (Random 
Two), and 1.03 (Random Three) times as good as the best conventional search strategy. The 
Random One persistence gf/hf ratio was 15/11; the genetic algonthm produced 
chromosome was f1e90234. The Random Two gf/hf = 15/4; GA-produced chromosome = 
f0b947c1. The Random Three gh/hf = 15/5; GA-produced chromosome = f1f73351. 
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Figure 24 : Random One Results 
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Figure 25 : Random Two Results 
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Figure 26: Random Three Results 
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2. Random Four 

This terrain was significantly harder than the other three randomly generated 
terrains which can be observed by the low performance of the search strategies. The 
average fitness of all strategies was 36% (Figure 24). The difficulty probably comes from 
the encapsulation of the goal. Examining figure 33, page 46, we can see that the goal is 
blocked by mostly high density blocks from (11,15) down to (11,10) across to (15,10). The 
only passible blocks are (15,10) and (11,12). Neither of which are a direct route, 
necessitating significant exploration. The genetic algorithm produced (f0b51535) scheme 
was 1.10 times better than the next best which was a persistence strategy with a gf/hf ratio 
of 14/15. This again seems to suggest that the genetic algorithm is only required when the 


problem is hard. 
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Figure 27 : Random Four Results 
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3. General Comments 


It is difficult to improve on the simple search strategies when the terrain is of low 
complexity. As suggested by DeJong [De92} the genetic algorithm can only optimize to a 
certain point (dependent on implementation) before reaching a state of dynamic 
equilibrium. The first three random terrains were of insufficient complexity to allow the 
genetic algorithm to convincingly surpass all conventional search schemes. It was however, 


in all cases, better than the best conventional ones. 


C. GENERAL COMMENTS 


In all cases, although the genetic algorithm produced strategy was always as least as 
good as the next best, it was not a substantial improvement over Persistence search unless 
the terrain was natural. Only in the most complex of the four random terrains did the genetic 
algorithm produced scheme really excel. This seems to suggest that the additional 
heuristics are only essential in natural terrains where some pattern in obstacle density exists 
or in random terrains of high complexity. 

Actual natural terrains, although usually best modeled by our natural terrains, could 
possibly be more similar to the random. Since the genetic algorithm produced search 
strategies are substantially better for our natural terrains and as least as good as standard 
search schemes for random terrains, they should be advantageous to use on any actual 
natural terrain. This is of course contingent on the physical agent’s dependence on minimal 
steps and its computational speed. If it’s computational speed is sufficient to avoid delays 
before each step and/or minimal steps are essential, the genetic algorithm produced scheme 
should always be used. 

Appendix C shows a comparison of the average time required for each strategy to 
search from start to goal for each of the terrains. As expected, the more complicated 
Strategies require additional computation time, but are not considered slow enough to 


prohibit their use except in cases of high speed agents with slow computational speed. 
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IX. CONCLUSIONS 


Heuristics previously used for search of an unknown space by a physical agent are 
distance from goal and distance from current. These are insufficient to minimize energy 
expenditure (steps taken) when some general knowledge of the area is known. The 
additional heuristics found to be pertinent are distance from start, crowding factors which 
account for obstacle node density around the considered frontier node, move-away factor 
which encourages reduction of the search space, and momentum which avoids wasted steps 
in course variations. These seven heuristics with their proper individual biases were found 
to be superior to standard search schemes. In this thesis we showed that genetic algorithms 
can be effectively used to develop optimal! heuristic biases that are adaptable to unknown 
search spaces if some general knowledge of the search space is available. Training done 
with randomly generated search spaces having common characteristics lead to robust 
search schemes which are, on the average, more fit than previously used strategies. 

We believe that this methodology of identifying all possible heuristics, fitting them 
into a binary representation, and applying genetics-based training is also applicable to a 
multitude of real-time search/optimization problems. Tests in other specific areas are 
needed to prove our conjecture. In addition, further research could be done in the 
application of more advanced genetic algorithms. Our results showed significant 
improvement using only basic genetics-based concepts, advanced techniques should 


continue to improve the effectiveness of resultant strategies. 
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APPENDIX A 
TERRAIN DEVELOPMENT FROM A SAMPLE DENSITY MATRIX 


Terrains are randomly produced using a density matrix as a guide. Figure 28 shows a density 
| matrix that was used to develope the terrain is shown in figure 29. This density matrix was not used 
for our analysis, but helps to make clear the relationship between the density matrix and the actual 
terrain. 

The density matrix is stored as a text file as shown in the figure. At each cycle for training or 
iteration for testing the density information is used to form a new terrain. Each hexadecimal number 
represents the desired density for a 4x4 area. The actual obstacle placement is random. Compare 
figures 28 and 29. The top left 4x4 area was filled in by checking if a random number (between 0 
and 15) is less than 4 at each node. This should on the average happen 4 out of 16 times making the 
obstacle count of each 4x4 area equal 4. The top left 4x4 is the average case with 4 out of the 16 
nodes being obstacles. 

The remaining 4x4 areas are filled out in a similar fashion. The start and goal nodes are chosen 


at random in the (2,2) and (13,13) areas as also demonstrated in figure 29. 


4444444444444444 
4444444444444444 
4444444444444444 
4444444444444444 
4444444444444 444 
4444444444444444 
4444444444444444 
4444444444444 444 
4444444444444444 
4444444444444444 
4444444444444444 
000000000000 0000 
C0102 03040506070 
O000000000000000 
S0g0ad0boOcOdDe OED 
CocoCoO00C00000CODG 


Figure 28 : Sample Density Matrix 
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Terrain has a density of 4 in the top 11/16 Goal 


BEEES EEBERERERE ESESE EREESEER SEES 2S SSA Bee | 
saon0 = Semeenen | i | ee -— ,eeneeneees aoe a 
al te + +. = 


SEES EEE EEE EEE EEE EEE EEE EEE 
le GEESE GEERRis Eecen i; oe es | | ttl os 
[| 888 Baanr || REE EERE Gee | 
0 tt — —- ra) an He tt CaRaE "..n. or 
[peel IE aera EEE Eee a 

a CoP 
L ERR eee eee 
ASC .cb Soc SEE eee 
S000 (0000000000000 00 0000S 0cee oe eee eee eee eee 





Start The lower 5/16 shows increasing density areas 0 to 15 
surrounded by areas of 0 density. 


Figure 29: One of Many Possible Resultant Terrains 
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APPENDIX B 


RANDOM TERRAIN DENSITY MATICES 


HMHWOAMOAINANHNMNMHM O 
HOM QWATTNQOEMW QM 
COQANM COM Ge AWDON ONM 
aAOMUDGEr OTe OM ~ornwrTim 
AAMHr-MROONTwWWdVO?T dE TG 
AMM OQAMUDBDHWHOTDTMNQWNAHM 
INON ST ST BGOrR DQ TC¢q¢ Nn HH 
HHOTATPAUNMNWOWTMWKWHM 
ANATOWOWM™ HONTIS © FB Bio 
WNDWDATIMNW QWONMNWMYT OY 
MHAMDDOOMTIwW &Y TWO ONTO 
HOMNOMdTM OW OTI1YKRK ¢ 
YOUWUTAGHBAMNOUGH INCOM 
AQAONMN DOT OHWMN ONAN 6M 
N COWD OM OT1 FTQH OWOTG 
GAT OHNMNNDIOMNNTNANMHN 


Figure 30: Random Terrain One 


MNVYVBUNOWMONRWnHNWMNwWSG 
YM TNMs CHR A gd HO 
WGTOAN THUMM CWOONnNWAMM 
Tomo QHnMnvssonn cnMmM 
ONTPTwOT AGH FsonTIMoMm 
BSENUNDWOANM QMOAME™ WY 
GUNMONHN DW YODA WrtMrico 
GCNOWVODVINH Grin CHD 
NTPOWPAMAWNNA FOM O 
COMD dD VYOWUNRMBM ww 
OMNMOANUDGTOAMNTIOSO GQan 
GNOHMMUBOMTOTOTON 
OAHWHTHr WO OWMNnHnn Aww 
WM QWnIMAaATtOoOM OWNDOT w& 
THOT1SC FT MUBUMNANMAWHM 
WtrtitWOAnMDMMNTOUVAT ITN 


Random Terrain Two 


Figure 31 : 
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SNINW O GHYMUNM sr BTS BMH 
NIN VDGNMWM OYMMW OWM 
YNDOQNMKNMNHGWMOAOMNOTD A 
COFTFTMHMRMHGDY OM Caer 
CWO SODVHNnDGGTSOVDGATM 
YOWUW AT ONWRHA CNOMY 
FANWHOMNANYE HR WO NnMNO 
YToO BODO GBMWONWdAA 
VAT VHP MINMIN DHNONY 
CODMDMDD BCDB ONDO VHRMNrewWM 
OMFrMWMNNS CHNYT ONM GG 
MMA ON QeuMnotvrsT Oo COOWw 
Ter yY OQVIODOO I MNnNMNroOdTCn 
TOM ODTAWS OMY INOMA 
CONQONTNOEXDHNAQWK YT ww 
NOTIM QNQTNNS Trim oOON 


Figure 32 : Random Terrain Three 


OBGDMUBMNNPFHO GTOAMNG QD G 
FANTIAMO QOMIBGWOM CMM 
TOWH GInig¢ CTHONYT TTL 
ANNE MW GHGiNnNon suHm 
AGSMNGslwods Ce Bio Grins 
Grn GMnrningoMmMIgwH7aA 
Ormoodyr onnwonaan 
NONE PTFWOMDY OYHnond eH 
HroHouumMmarnHodgtag @ 
NWOOMT1 GCNHOOoT GH ow 
ANTINWYUHUMRNMOODMDWOOOCO Ar 
OVUM OMNWDWHTFO OMNMNTS OM 
NOMnNGT CNAMYDTMODOOWW 
IN BCOWWD YOU DTM Gonw os 
YOO OQronwiwnygyaAnwdoon 
AMNWEMNMTINOAMoODG IM 


Figure 33 : Random Terrain Four 
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APPENDIX C 


SEARCH STRATEGY COMPUTATIONAL TIME 


TABLE 1: AVERAGE SECONDS REQUIRED TO SEARCH EACH TERRAIN 


; ; Real-Time- GA- 
Central 
Mountain 0.0239 0.0274 0.0065 0.0163 0.0615 
Single Left 
Ridge 0.0232 0.0264 0.0075 0.0176 0.0694 
Singe Right 
Ridge 1.3051 0.1496 0.0199 0.1405 0.2718 
Double 
Ridge S253 1.3562 0.0420 0.4131 2.3420 
| Single Left 
Plateau 0.1071 0.0982 0.0167 0.0626 0.1481 
Single Left 
0.0769 0.0856 0.0146 0.0487 0.1310 


Plateau 
With  _— 
Random 
One 0.0177 0.0265 0.0114 0.0478 
Random 
Two 0.0313 0.0267 0.0056 0.0150 0.0644 
Random 
Three 0.0294 0.0258 0.0065 0.0151 0.0521 
Random 
Four 0.1629 0.2255 0.0165 0.048 1 0.1563 
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le: 


Zee 


APPENDIX D 


PROGRAM C CODE 


Ga (seater e @#*eoeoeee#egoe ® e e@ e@ e@ e@ e e @ e 

mala .c¢ e co e e oe e e e ee eeeoeee ® e ® e e e e e e 
main 

trainm.Cc.... =e 7 ae 
train 
put_gen 
put_rs 

tesc..c e@ e e e e@ oe e@®eoeee#e#e#e#eeeee® e e e e e e 
test 

ESCEUP ae e e@oeoeeeeeoeneee#eee ee @ e e e e@ e e @ 
gets seed 


get persrch-om 
read_density_ file 
make array 

make node 

find node 


tpopu Ta.b..On eC a ee ee 


create population 
new_individual 
get population 
put population 


astar.Cie... .2 eee. aoe 


a star 
update_astar frontier 


oie esr lere Sp oh Go SRI 6 
ete Clams 
move_adjacent 
find best 


YC gas Cale, C oe ee i ote alec eee 


ei ae aie i 
update _rtastar_adjacent 
insert 


bisearchic ..2.0gee Pe Me sna 


bfsearch 
bf update frontier list 
Bf pick bestatronticr 


pseaween.¢ see, ee ee 


psearch 
p_update frontier list 
p_pick_best_frontier 


CSEaArCN.cC .. 52 Pr ee ee eae x 


search 
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m3 . fominIa@Gii Ol . C4... cs sees <s6 6 8 eS evsce © oe 4060 CeCe 4 gow ee ele ee « a4 


update frontier list 
update adjacent_obstacles 
pick _best_ frontier 
reset_back_ track_state 
update_list 

ehigers plays 

adjacent 

update_crowd_ sides 
update crowd diag 
calc_move_away 
calc_momentum 
compute_subtotal 


14. CLES SNC. pie 5 Pa ae Gh ae ©. & Se 8 4 


insert_heap 
delete heap 
move_heap 
swap 


eS . SyOmVve.C .. ks ees ~ GaSe ee Separates sine’ Queene 104 


evolve 
eredaecemacenneap 
gamete Crossover 
bit crossover 
crossover 
get_mask 

mutate 

one_if mutate 
get_odd 
set-equal 


EO SENS Ch es a ee eee eee ee es 


insert _mate_heap 
pop_mate_ heap 
move_mate_heap 
Swap num 





ey . Cn Ge 5 eR 5 eeeee?#e # @# @# @® @® @ Sly. in ee de 4 
gen_xi 
gen_yi 
equalf 
show_least_nodes 


1s. rai goals, © a ee ea wei eel 7 

initialize 

draw_terrain 

show_mouse 

draw _nodes 

draw grid 

squaref 

square 


a9 , (eC. ok 5 ee eee#e#ueeoeo#e#»n°e?eeteeeeee ee # @ ee Ses cee LL 
print_density 
print_population 
print_node 
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ga_search.h 


/* 
File: ga_search.h 
Programmer: g.b. parker 
Environment: any 


Language: & 

Date: 9 july 92 

Revised: 

Comments: This file holds all the header information needed for all 
ga_search files. 


a7 
#include <stdio.h> 


#define 


0 
#define 1 


Ee 
A 
#define UNTOUCHED 
#define OBSTACLE 
#define VISITED 
#define FRONTIER 
#define START 
#define GOAL 
#define CURRENT 


#define SHORTEST 
#define X 


Or DO B&B WN FE © 


#define NUM 0 


#define N 
#define E 
#define § 
#define W 


#define NE 2: 
#define SE 2 
#define SW 3) 
#define NW 4 


#define MASKO 0 /* 0000 */ 
#define MASK1 8 J* 10007 -*7 
#define MASK2 12 JX TAO Ox 
#define MASK3 14 J®-LY1O 7, 
#define MASK4 15 /* (Adee, 
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#define NA =1 /* not applicable; for frontier in- 
dex */ 





#define BIG NUMBER 10000 
#define STANDARD DENSITY 4 
#define PROB BIT MUTATE 50 /* x/10000 prob of mutate */ 
#define SQRT2 1.414213562 
#define randl6() ((random()/13) % 16) /* return rand int from 0 to 15 */ 
#define rand8() ((random()/13) % 8) 7* return randwanteGromm0. © Ome, 
#define rand5S() ((random()/13) % 5) /* return rand int from 0 to 4 */ 
#define rand1l0000() (random() % 10000) /* return rand int from 0 to 9999 */ 
/* Graphics definitions */ 
#define SLEEPTIME 2 
#define ASTITLE vA Search 3 
#define GATITLE "Genetic Algorithm Produced Search" 
#define PERSTITLE “Persistence Search " 
#define BFTITLE "Best First Search 
#define HCTITLE Mati Climb Search : 
#define RTASTITLE “Real Time A* Search iy 
/* Node record for terrain */ 
struct node rec 
{ 
nt Xi; 
nit Yi; 
EvOat xX; /* for graphics */ 
mloat y; /* for graphics */ 
int state; /* UNTOUCHED, OBSTACLE, VISITED, or FRONTIER */ 
amt back track state; /* UNTOUCHED, VISITED, or OBSTACLE */ 
float subtotal; /* includes all but dist from current & move away */ 
Bloat Gist from start; 
Proagt Gist #from¥goal; 
float dist _from_current; 
struct node rec *predecessor; /* points to predecessor for a-star search */ 
int frontier index; /* position in frontier heape*/ 
Struct node_rec *qnext; /* next in q for dist fromecurrent DFS */ 
struct node _rec *qreset; /* reset link list after back DFS */ 
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} 
/* Individual record for population */ 


Struct factor Sseruccay 
unsigned int place holder 
unsigned int start dist 
unsigned int goal _ dist 
unsigned int current_dist 
unsigned int crowd_sides 
unsigned int crowd diag 
unsigned int move_away 
unsigned int momentum 


aie 


e =e ‘“e 


‘e ™e ee 


a» bf bh bf fb Lh DL LHL 


‘ae 


union chrom_union { 
SUEWEE [actor Struce staceon: 
unsigned int alleles; 

boG 

Struct andivadual eseructe.| 
union chrom_union chrom; 
int fitness; 
eAKeyele satis Ch Cs 
int previous index; 


ae: 


/* global variables */ 
extern int heap _ size; 


/* Functions listed under file */ 


{/* astar,c */ 
float a star (jy 
int update _astar frontier(); 


jemeemastar.c */ 

Proatrrta star (), 

int update_rtastar_adjacent(); 
int insert(); 


/* hill seieimp Fern / 

Float isihine emo, 

struct node_rec *move_adjacent(); 
Struct node_rec *find best (); 


/* test.c */ 
int test (); 
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P= train.c */ 
ant train(); 


¥* tsearch.c */ 
float search(); 


/* psearch.c */ 

float psearch(); 

ait ep update frontier list (); 

Meret enoce rec *p pick_best._frontier(); 


/* bfsearch.c */ 

float bfsearch(); 

mit bf update frontier list(); 

struct node_rec *bf pick best _frontier(); 


/* tsetup.c */ 

int get_seed(); 

unsigned int get_pers_chrom(); 
int read density file(); 

int make _array(); 

int make _node(); 

Struct node rec *find node(); 


/* tpopulation */ 

int create population(); 

struct individual struct *new_individual(); 
mt get population (); 


fe Eprint.c */ 

mt print density (); 
int print _node(); 
Piteprant population (); 


emecrcontier.c */ 

int update frontier list(); 

Seruet: node rec *pick best_frontier(); 
int update _adjacent_obstacles(); 

float update dist (); 

int update crowd sides(); 

int update _crowd_diag(); 

int calc_move_away(); 

int calc_momentum(); 

float compute _subtotal(); 
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/* theap.c */ 

int insert heap(); 
int delete heap(); 
int move_heap(); 
int swap(); 


/ *emiscues +7 

float compute shortest); 
int equalf(); 

/* update dist stant (79 */ 
int gen_xi(); 

int gen_yi(); 

int show_least_nodes(); 


/* evolve.c */ 

int evolve(); 

int create _mate_heap(); 
int crossover(); 

int allele crossover (); 
Int “Bit SCLOssover ()i; 
int get _mask(); 


/* eheap.c */ 

int insert _mate_heap(); 
int pop mate _heap(); 
int move_mate_heap(); 
int swap _num(); 


/*tdisplay.c */ 

int initialize(); 
int draw terrain(); 
int show_mouse(); 
int draw _nodes(); 
int draw _grid(); 
void squaref(); 
void square(); 
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| 


main.c 


/* 
File: main.c 
Programmer: g.b. parker 
Environment: any 


Language: e 

Date: 9 july 92 

Revised: 

Comments: titsis the Control for user input and call of train or test. 
Ten command live arguments are optional. The syntax for a call is as follows 
(the 0 argument is the program call): COdiss Cee ()) OLfor «rain. |). for tese, 


(2) random seed; (3) input population file name; (4) input terrain density file 
mame; (5S) start region on X axis; (6) start region on Y axis; (7) goal region 
on X axis; (8) goal region on Y axis; (9) number of generations if training, 
iterations if testing; (10) cycles per generation if training, array position 
of best individual in the GA_produced population for testing; (11) file name 
for out population if training, hexadecimal representation of best Persistence 
search scheme. 


e/ 


#include “ga_search.h" 


/* KKK KKK K KKK KKK KEKK main KKK KKK KKK KKK KKK KK KKK KKK KEKE KKK KEKE KKEKKKEKKKKKKK x / 


‘main( argc, argv ) 


mit argc; 
Char *argv{]; 
{ 
Beeoct individual struct *individual (32); 
int arg _ seed = 0; 
euar arg Population[32]; 
char arg population out [32]; 
char arg density([32]; 
int choice = 0; 
ie sx = 2; 


mie. sy = 2; 

int gx = 4; 

ie ogy = 2; 

ime iterations = 25; 

int generations = 3; 

int cycles per generation = 2; 
int best_individual = -1; 


unsigned int pers chrom; 
char hname[64]; 
int hnlength; 
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gethostname( hname, hnlength ); 


Strcpy( arg population, 


strcpy( arg population_out, “popx.out 


strcpy( arg density, 


Switch ( argc ) 
{ 

case 12: 

sscanf( argv[l1l], 

sscanf( argv[ll], 
case ll: 

sscanf( argv[10], 

sscanf( argv[10], 
case 10: 

sscanf( argv[9], 

sscanf( argv[9], 
case 9: 

sscanf( argv[8], 
case 8: 

sscanfi( argv(?])), 
case 7: 

sscanf( argv[6], 
case 6: 

sscanf( argv[5], 
case 5; 

sscanf( argv[4], 
case 4: 

sscanf( argv[3], 
case 3: 

sscanf( argv[2], 


be fs fae 
Ww tx" 


we =a. . 
8 Eel 


wx Gl - 


Salo 
aie (ol 
ie fel te 
a oclae, 
"tae, 
"3s", 
eS 


we qa" ; 


arg populations oue.).; 
&pers_ chrom ); 


&cycles per generation ); 
&best_individual ); 


&generations ); 
&iterations ); 


&gy ); 
&gGx ); 
&sy ); 
BS), 


arg density ); 
arg population ); 


&arg_seed ); 
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ey 
MG 
Ves 


case 2: 
sscanf( argv{1), "%d", &choice ); 
switch ( choice ) 
{ 


case 0: 


train(individual,arg_ seed,arg_population,arg density, sx, 8y,gx,9y, 


generations, cycles_per_generation,arg population_out,hname ); 
put_population( individual, arg population _out ); 
break; 


case 1: 
individual[0Q]) = NULL; 
test (individual,arg seed,arg population,arg density, sx, s8y,gx,9y, 


iterations, best_individual,pers_chrom) ; 
break; 


case 2: 


train(individual,arg seed,arg population,arg_ density, sx, sy, 9x,9Qy, 
generations, cycles per generation, arg population_out, hname) ; 

test (individual,arg_seed,arg_population,arg density, 8x, 8y,9x,9y, 
iterations,best_individual,pers_ chrom) ; 


put _population( individual, arg population _out ); 
break; 


} 
break; 
case 1: 
train(individual,arg_ seed,arg_population,arg density, 8x, S8y,9x,9y, 
generations,cycles per generation,arg population_out,hname) ; 


put_population( individual, arg_population_out ); 
} 


oY 


train.c 


/* 
File: tha liee 
Programmer: g.b. parker 
Environment: any 


Language: S 

Date: 9 july 92 

Revised: 

Comments: Called by main to train a population of 32 individuals.) )teege 
input population, a random one is generated. 
*/ 


#include "ga_search.h" 


/* KKK KKK KK KKK KK train KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK K x / 





train( individual, arg seed, arg population, arg density, sx, sy, 9x, gy, gen- 
erations, cycles per generation, arg population_out, hname ) 
struct individual Struct “indir dwar 22, 
int arg_seed; 
char arg population([32]; 
char arg _density[32]; 
int. Sx, SY ,0gx, aos 
int generations, cycles per generation; 
char arg populattenmeur [321 
char hname([64]; 
{ 
int dens rey (tea vo); 
StrUCE Node r rec. “Node (co lmloo iF 
int gen, cycle, a, rs, Short scoune, 
Float shortest path; 
int dummy; 


rs = arg seed ? arg_seed : get_seed(); 
srandom(rs); /* seed the random generator */ 


printf ("\nRandom seed is %d", rs ); 
get_population( individual, arg population ); 


read density file( density, arg density ); 
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for( gen = 1; gen <= generations; gen++ ) { 
Brinte(™\n gen = %*d (cycle,rs) — “, gen 
for( cycle = 1; cycle <= cycles per generation; cyclet++ ) { 
rs = rs + 1; 
Sone count = 0; 
while( (shortest path = a_star(sx,sy,gx,gy,rs,density,node )) 
> ( BIG NUMBER - 1.0) ) { 
if( short count > 1000000000 ) { 
printf ("\nPROGRAM ABORTED - iteration %*d - no shortest path\n",i); 
ECturN{(F); 
} 
else 
Poe— 7 To + 1; 
} 
Preiget(" “(g0,td)", cycle, rs ); 
for( 1=0; 1<32; i++ ) { 
if( cycle == 1 ) 
individual([i])->fit_sum = shortest path / search( sx,sy,9x,9y, 
individual [i]->chrom. factor, rs,density,node, &dummy ); 
else 
individual[iJ->fit sum = individual[i])->fit sum + shortest path / 
search( sx,s8y,g9x,gy, individual [i]->chrom.factor,rs,density, 
node, &dummy ); 
if( cycle == cycles per generation ) 
icivedual|ij—=>fitness = (int) ((individual[i]->fit_sum / 
cycles per generation) * 100.0); 


} 
if (gen == generations - 1 ) 
cycles per generation = cycles per generation + 10; 


evolve( individual, rs ); 
puewrs( rs ); 
if ( ( (gen % 10) == 0) II (gen < 10) ) 

put _gen(gen,arg population _out,hname, rs) ; 

/* put gen to a standard update file */ 

if ( (gen % 50) == 0) { 

put population( individual, arg population_out ); 
} 
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/* KKK KKK KKK KKK KKK put gen KKK KKK KEK KKK KKK KKK KKK KKK KKK KKK KKK KKK x / 


/* Called by train to continually store status information to a file in the 
directory of execution */ 


put_gen( gen, arg_population_out, hname, rs ) 
int gen; 

char gq populatwon coum); 

char hname[64]; 

Iintess, 


{ 
FILE *gen_ file, *fopen(); 


gen_ file = fopen("running.update", "a"); 


fprintf(gen_file," %s *s gen = %d rs = %d\n", hname, arg population out, 
gen, rs); 


fclose(gen_file) ; 


/* KaK KKK KEKKKKKKKK KK put rs KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKKEKKK x / 


/* Puts random seed info to a file in the directory of execution */ 
put. rs( ose) 
Lint Es, 
{ 
FILE *pegtile,), *fopeniyy 
rs_ file = fopen("rs.update", "a"); 


mopanbevens (ity je atiley Ue sele bee})) 6 


felose (rs tiie); 
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test.c 


/* 
File: test.c 
Programmer: g.b. parker 
Environment: any 


Language: : 

Date: 9 july 92 

Revised: 

Comments: Called by main to perform a comparative test of search 


Strategies. The default is for all individuals of the population to be 
tested, unless a specific individual is specified. 


e/ 


#include “ga_search.h" 
#include <sys/time.h> 


/* Kea KKK KKK KKK KKK K test KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KK * / 


test ( individual, arg _seed, arg population, arg _density, sx, Sy, 9X, Jy, 
iterations, best _individual, pers_chrom ) 

mercer individual struct *individual [32]; 

int arg _ seed; 

Ear arg population([32]; 

char arg density[32]; 

mit sx, SY, GX, 9y; 

mit iterations; 

int best_individual; 

unsigned int pers_chrom; 

{ 


Seruct individual struct *best first; 
PeEUcceindividual Struct *peérsistence slearch, 
int density[16] [16]; 

Struct node rec *node[66] [66]; 


Mite, ©S, Kk, Short count; 


float shortest _path; 

float realtime astar_fit_sum; 
int realtime_astar_ fitness; 
wijdtamilleclimoa £1t sum; 

igen eelimb> fitness; 
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float temp; 

float gat = 0.0; 
float gavrcks = 
PCat Omit =U), 


float bE ticks = 020, 
float pers t = 0.0; 
float pers ticks = 0.0; 


float Ne st. = 3080); 
float he e2#eks = 9070, 
float rea t = 020; 
float rtageveks= Geer 


long sec, usec; 
struct timeval *tvp = (struct timeval *)malloc(sizeof(struct timeval)); 
struct timezone *tzp = (struct timezone *)malloc (sizeof (struct timezone) ); 





best first = 
(struct individual struct *)malloc(sizeof (struct individualiscteuer eae 
persistence search = 
(struct individual struct *)malloc(sizeof (struct individual )seeue me 


best _first->chrom.alleles = 0xe0100000; 
persistence search->chrom.alleles = pers chrom; 
best first—>£it (sum = 0.0; 


persistence search->fit sum = 0.0; 


ES) == aro SCC esd ge see Germs e—c nr, 
srandom(rs); /* seed the random generator */ 


printf("\nRandom seed is %d", rs); 
if ( individual[0] == NULL ) 
get_population( individual, arg_population ); 


for (k=0;k<32;k++t) 
individual [k]->fit_ sum = 0.0; 


read density file( density, arg_density ); 


PriINtt (7 \n (iteration, Ss) oe 
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Bort oa = 1; 1 <= iterations; i++) { 
Pio s + 1; 
Short count = 0; 
while ( (shortest_path = a_star( sx, sy, gx, gy, rs, density, node )) 
> ( BIG NUMBER - 1.0) ) { 
taieshort count > 100GQ000000 ) { 
printf("\n PROGRAM ABORTED - iteration %d - no shortest path\n", i ); 


Return (rl); 

} 

else { 
Pemer™ (td,sdy",i,rs ); 
rs = rs + 1; 


} 
Peatmen( (sd, sd) ",1, rs.) ; 
if (best individual == -1) 
for (k=0;k<32;k++) { 
temp = search( sx,sy,9x, gy, individual [k]) ->chrom. factor, rs,density, 
node,&ga_t ); 
individual [(k]->fit sum = 
individual [k]->fit_sum + shortest path / temp; 
if (temp < shortest path ) 
printf ("\nSHORTEST PATH > ACTUAL PATH "); 
} 
else { 
k = best_individual; 
temp = search( sx,sy,gx,gy, individual [(k]->chrom.factor,rs,density, 
node,&ga_ t ); 
individual (k]->fit_sum = individual(k]->fit_sum + shortest_path / temp; 
if (temp < shortest path ) 
Printe¢ \nNSHORTEST PATH > ACTUAL PATH "™); 
Fcgselekse= ga ticks + gait; 
best first=->fit_ sum = best first->fit_sum +. shortest_path / 
bfsearch( sx, sy,gx,gy,best_first->chrom.factor, 
rs,density,node,&bf t ); 
Brewicks = bE ticks + bf t; 
persistence _search->fit_sum = persistence_search->fit_sum + 
shortest path / psearch( sx, s8y,9x,gy, 
persistence search->chrom.factor, rs,density,node,&pers t ); 
pers ticks = pers ticks + pers t; 
hill climb fit sum = hill climb fit sum + shortest_path / 
Pec linntmsx,—S¥,eOx,9y,6Ge, density, node, &hc t ); 
emt ekss=ehemercks + he t; 
realtime_astar fit _sum = realtime_astar_fit_sum + shortest_path / 
rta_star( Sx, sy, gx, gy, rs, density, node, &rta_t ); 
CedmerGkoemmrta ticks + rta t; 
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if (best_individual == -1) /* ne best andhvadia Seer ce, 
for (k=0;k<32;k++) 
individual [{k]->fitness = 
(int) ((individual[k]->fit_sum / iterations) * 100.0); 
else 
individual{k]->fitness = 
(int) ((individual[k]->fit_sum / iterations) * 100.0)% 


best_first->fitness = (int) ((best_first->fit_sum / iterations) * 100.0); 
persistence_search->fitness = 
(int) ((persistence_search->fit_sum / iterations) * 100.0); 
realtime_astar fitness = 
(int) ((realtime astar fit sum / iterations) =~ 100m0ne 
hill climb fitness = (int) ((hill_climb fit sum / iterations) * 100.0); 


if (best individual == -1) 
print population( individual ); 
else { 


printf("\n ga-produced %2d %7.3f %7.3f£ $x", individual[{k]->fitness, 
individual [k]->fit_sum, ga_ticks7 individual [ki=>chrom-alleles @ 
Printf£ ("\n best sfixrst ——-s2de27 38 ewer tx", best _first->fitness, 
best _first->fit_sum, bf _ ticks, best_first->chrom.allteles); 
printf ("\n persistence %2d %7.3f£ %7.3f Sx", 
persistence _search->fitness, persistence_search->fit_sum, 
pers ticks, persistence search chron jd lleres:. 
printf ("\n hill climb 2d %7.3£ 47 .3f", hieeePimneeeene ss 
hill clamb f£2tusum, he scickses, 
printf("\n RTA star S20 47.3f %7.3£", realtime astar fitness, 
realtime astar fit sum, rta_ticks ); 
} 
Printt( Ale: 
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tsetup.c 


/* 
File: ESecLUup.c 
Programmer: g.b. parker 
Environment: any 


Language: e 

Date: 9 july 92 

Revised: 

Comments: Setup functions 
e/ 


#include "ga_search.h" 


/* kkk kkk KKK KK KKK KKK get seed KK KKK KKK KK KKK KKK KK KKK KKK KKK KK KK KKK KKK x 


/* Interfaces with user to get random seed */ 


get_seed() 
{ 
char nl{l]j: /* absorbs new_line after seed entry */ 
int rand_seed; 
printf("\nEnter random seed or 0 (system assign seed): "); 
scanf ("td", &rand_ seed); 
gets (nl); 


if (rand_seed != 0) 
return rand seed; 
else 
return getpid(); 


/* kek kkk keke KKK KKK KKK get pers chrom waa KKK KKK KAKA KKK KKK KKK KK KKK KKK KKK K x / 


/* Interfaces with user to get Persistence chromosome */ 


unsigned int get_pers_ chrom() 


{ 


Ghar nij{l)j:; /* absorbs new_line after seed entry */ 
unsigned int pers chrom; 

printf("\nEnter Persistence chromosome (8 hex digits) or 0 
scanf("%x",&pers chrom) ; 


gets (nl); 

if (pers_chrom != 0) 
return( pers chrom ); 

else 


return( 0xe0110000 ); 
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(e0110000): 


as 


[& kk KKK KKK KKK KKKK read density File *%%—KK III II III I / 


/* Reads density file from execution directory */ 


read density file( density, file_name ) 
int density[16] [16]; 

char file name[32]; 

{ 


nih 9 ei aRp eam 6: 

FILE *density file, *fopen(); 

int not_end = T; 

int node_density = STANDARD DENSITY; 


if( file _name[0] == ' ') { 
printf("\nEnter density file name: "); 
gets(file name); 


} 


if ( (density file = fopen(file_name,"r")) == NULL ) { 
printf("\nThe file does not exist, standard densities being used."); 
not end = F; 

} 

for ( j=15 ; j>=0 ; j-- ) { 
for (-i1=80 (7 gisele ee 


if ( not_end && ( fscanf(density file,"%x",&node density) != EOF ) ) 
density[iJ][j] = node density; 
ele 


density[i][j] = STANDARD DENSITY; 
not_end = F; 


} 
} 


if (density file != NULL) 
fclose(density file); 
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/* KKK KKK KKKKKKEKKKKK make array 


KOK KKK KKK KKK KK KKK KKK KKK KKK KK KKK KKK KKK x / 


/* Creates the node array on initial use, 


make _array( density, node ) 
eat density(16](16]; 

struct node rec *node[66] [66]; 
{ 


static int first T: 


mt 1, Vs 

for (i=0;1i<=65;it+) { 
make node( node, 
make node( node, 
node[(i](0]->state 
node[i] (65]->state 


0, 
65, 


i, 
ly 


One ls j<=64;5++) { 
make node( node, 
make node( node, 
node [(0](j3]->state 
node [65] [(j3]->state 


0, 
65, 


j, 
J, 


for (j=1; 3<=64; j++) 

for (i=1;i<=64;it+) { 
make node( node, 

if (randl6() 


i, 


node[{i}][{j]->state 


J, 
< density[(i-1)/4] ((3-1) /4]) 
OBSTACLE; 


fer Ste 


f.rse); 


OBSTACLE; 
OBSTACLE; 


first); 


fxrst )e: 


OBSTACLE; 
OBSTACLE; 


Er Sty )y: 
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then resets records after that */ 


/* indicates if first time to make array */ 


/* KKK KKKKKKKKKKKEK make node KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKEKKKKK x / 


/* Resets single node information */ 


make node( node, xi, yi, first ) 
struct node_rec *node[66] [66]; 
Int. so 

ngle, 33731, ¢ 

int first, ee) * eer Pe eX 


{ 


int k: 


npG (a oealietce, )) weal! 
node[xi] [yi] = (struct node_rec *)malloc(sizeof (struct node_rec)); 
node [xi] [yi]->xi = xi; 
node [xi] [yi]->yi = yi; 
node [xi] [yiJ]->x = (float) xi; 


node [xi] [yi] ->y = (float) yi; 
} 
node [xi] [yi]->state = UNTOUCHED; /* OBSTACLE, VISITED, or FRONTIER */ 
node [x1] [yi] ->back_track state = UNTOUCHED; /* VISITED, or FRONTIER */ 


node [xi] [yi] ->subtotal = BIG NUMBER; 

node [x1] [yi]->dist_from_start = BIG NUMBER; 
node [xi] [yi] ->dist_from_goal = BIG NUMBER; 
node [xi] [yi]->dist_from_current = 0.0; 

node [xi] [yi] ->predecessor = node[xi] [yi]; 

/* points to predecessor for a-star search */ 
node [xi] [yi]->frontier index = NA; /* not have index to frontier heap */ 
node [xi] [yi] ->qnext = NULL; 

/* points to next in q for dist from current DFS */; 
node [xi] [yi] ->qreset = NULL; /*x reset link list after back DFS */ 
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/* KKK KKKKKKKEKK find node KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKK x / 


/* Picks a random node in the designated density area. Used to identify 
Start and goal nodes. */ 


struct node_rec *find_node( node, dens_col, dens_row) 
Struct node rec *node[66] [66]; 
int dens_col; 
_ int dens_row; 
{ 
int k, xi, yi, base_x, base_y; 


base x = (dens_col * 4) + 1; 
base y = (dens_row * 4) + 1; 


for (k=0;k<100;k++) { 
xi = base x + (randl6() % 4); 
yi = base_y + (randl6() % 4); 


if (node[xi] [yi]->state == UNTOUCHED) 
: return node [xi] [yi]; 


} 
return node([xi] [yi]; 
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tpopulation.c 


/* 
File: tpopulation 
Programmer: g.b. parker 
Environment: any 


Language: &: 

Date: 9 july 92 

Revised: 

Comments: Functions dealing with population creation/storage 
a 


#include “ga_search.h" 


/* KKKKK KKK KKEKKKKKKE EK create population KKK KK KK KKH KK KKK KEK KKK KKK KKKKEK x / 


/* Generates a population of random individuals */ 


create population( individual ) 
struct individualist ruct. *inerviduaeioe ie 
{ 

Sele | 


for (k=0;k<32;k++) { 
individual[{k] = 
(struct individual _ struct *)malloc(sizeof (struct individual struec) me 
new_individual( individual[k] ); 





} 


/* KKK KKK KK KKKKKKEKK new individual KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKEK x / 


/* Sets initial values of individual records fields */ 


struct individual struct *newsrndivicnalijjenas, 
Struct) mndividva lest rict a. tue. 
{ 

ind->chrom.factor.place holder = Oxf; 


ind=>chrom. factor. stant serait = randl6(); 
ind=>cnrom. faceor doa leas st = randl6(); 
ind->chrom: factor .current idist.—seanmclogs, 
ind->chrom.factor.crowd sides = rand1l6(); 
ind=>chrom.factor.crowdidiag §)— canara, 
ind->chrom.factor.move_ away = randl6(); 
ind->chrom. factor.momentum = randl6(); 
ind->fitness =O. 

ind=> fit esum = 8.05 
ind->previous_index = 99; 
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[& KKKKKKKKKKKKKKKKK Get population 
/* Reads population froma file */ 


KEK KKKKEKKKKKKKKKa KKK KKK KKK KKK KKK x / 


get_population( individual, file name ) 
Beruct individual struct *individual [32]; 
char file name[32]; 


{ 


Bho s our 

FILE *population_ file, *fopen(); 
Pitenot end = T; 

unsigned int alleles; 


if ( file name(0) == ' ') { 


printf("\nEnter population file name: "); 


gets(file name) ; 


if ( (population_file = fopen(file name,"r")) == NULL ) { 
printf("\nThe file does not exist, random population being used."); 


create population( individual ); 


} 
else { 


morte) 7 i1<32 7 it+ ) { 
moarvaadual[i] = 
(struct individual struct 
individual[i]->fitness 
mi@iviaual[i)=->f£it sum = 
individual[i]->previous_index 


if ( not_end && ( fscanf (popul 
individual[i})->chrom.alleles 


else { 


eyialloe(sizeotistruct individual struct) ) ; 


= 0; 


02a; 
= 99; 


ation file,"%x",&alleles) != EOF ) ) 
= alleles; 


new_individual( individual[i] ); 


not eena = Fb; 
} 


} 
fclose (population file); 
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* 
[ KRKKKKEKKKEKEKERK put Population  * ke RRR KR KEKE KEK KK KEKE KK KKK 


/* Puts the population to a designated file */ 


put_population( individual, file name ) 
Struct individualist ruct, *endividualmpe czar, 
char file name[32]; 


{ 


TG aane 
FILE *populatzon f1le, <Lopent): 
if( file name[{0]) == ° ' ) { 
printf ("\nEnter output population file name: "); 


gets(file name); 


} 
population file = fopen(file name,"w") ; 


for ( i=0 ; 1<32 ; i++ ) 
fprintf (population file, "%x\n", individual [i]->chrom.alleles) ; 


fclose (population file); 
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ai 


astar.c 


/* 
File: agiear.c 
Programmer: g.b. parker 
Environment: any 


Language: Cc 

Date: oF uly 92 

Revised: 

Comments: A-Star search - Finds the shortest path 
e/ 


#include “ga_search.h" 


int heap _ size; 


/* KEKE KKKKKKEKKKKKKK OB star KEK KEKE KEKE KEKE KKK KKK KKK KKK KKK KKK KKK KKK KEK KKK KKK x / 


float a_star( sx, Sy, gx, gy, random_seed, density, node ) 
ine SXx,SyY,9X,OY? /* position in density array for start & goal */ 
int random_seed; 
ent density(16] (16); 
Struct node _ rec *node[66] [66]; 
{ 
struct node_rec *current, *start, *goal;_ 
struct node_rec *frontier heap[4096]; 
ant kk; 


#ifdef IRIS 
union chrom_union dummy_cu; 
dummy cu.alleles = 0; 
#endift 


heap size = 0; 
srandom(random_ seed); /* seed the random generator */ 
make _array(density, node); 
Start = find_node(node, sx, sy); 
Sstart->state = START; 
Seatt—-dist from start = 0.0; 
goal = find_node(node, gx, gy); 
goal->state = GOAL; 
for (k=0;k<4096;k++) 

frontier _heap[k] = NULL; 
SUurrent. = Start; 
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#ifdef IRIS 
/* initialize the IRIS system */ 
initialize (ASTITLE) ; 
draw_terrain (node, start,goal,current,goal->dist_ from start, 
dummy_cu.factor); 


#tendif 


while ( current != goal) { 
update_astar frontier( node,current, frontier_heap,start,goal ); 


if( heap size == 0) break; 
current = frontier heap[0]; 
delete _heap( frontier heap, frontier heap[0} ); 


#ifdef IRIS 
draw_terrain (node, start, goal, current, goal->dist from start, 


dummy cu.factor); 
sleep( SLEEPTIME ); 
#endif 


goal->state = GOAL; 
Start->state = START; 


return( goal—-dist ftromistarea., 
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[®t kk IO update astar frontier  * ® % KKK RRR ITOK ITOK KOK IKK 


Medate astar frontier( node, c, frontier heap, start, goal ) 


SCrucct 
struct 
Struct 
struct 
struct 


{ 


oat Xl, 


base xi = c->xi 
base yi 
map Xi 

Repay. = C—>yi 


node rec 
node rec 
node rec 
Hoge rec 
node_rec 


cm>yi 
C->xi 


*node [66] [66]; 

“Cc; /* current */ 
*frontier heap[4096]; 
AStaxrt; 

Zaqod |; 


Wieoase x1l,—5Dase Yi, topsxia top yi; 


ite, 91 27 C=>Xing 1. 
eee CV — 
64 ? 64 : c->xi + 1; 
64 ? 64 : c->yi +1 


£Or (xi=mbase_ xi; xi<=top_xi;xitt) 
for (yi=base_yi;yi<=top yi; yi+tt) 


else if 


lp dee | 


(node [xi] [yi]->state == UNTOUCHED) | | 


(node [xi] [yi]->state == GOAL) ) { 
node [xi] [yi]->state = FRONTIER; 
node [xi] [yi] ->predecessor = ¢c; 


node [xi] [yi]->dist_from_goal 
node [xi] [yi]->dist_from_start 


update dist (node[xi] [yi],goal); 


Ga-arsenrtrom start + update dist (node[x1i] [yi],¢c); 


node [xi] [yi]->subtotal = node[xi] [yi]->dist_from_goal + 


node [xi] [yi]->dist_from_start; 


insert _heap( frontier heap, node[xi] [yi] ); 


( node[xi] [yi]->state == FRONTIER ) { 


if ( node[xi] (yi]->dist _from_start > 
(c->dist_from_start + update dist (node[xi][yi],c)) ) { 
node [xi] [yi]->dist_from_start = 


c->dist from_start + update dist (node[xi] [yi],c); 


node [xi] [yi]->subtotal = node[xi] [yi]->dist_from_goal + 


node [xi] [yi] ->dist_from_start; 


move _heap( frontier heap, node[xi]J [yi]->frontier index ); 


} 
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me) 


hill_climb.c 


/* 
File: hale Pimb-< 
Programmer: g.b. parker 
Environment: any 


Language: . 

Date: july. 

Revised: 

Comments: Hill climb search - Finds best adjacent frontier node or back- 
tracks 
4) 


#include "ga_search.h" 
#include <sys/time.h> 


/* KKKKKKKKK KKK KKK KK hill climb KKK KKH KKK KKK KKK KKK KK KKK KKK KKK KK KKK KKK x / 
float hill_climb( sx, sy, gx, gy, random_seed, density, node, ticks ) 
int. sx}sy,40*, ov /* position in density array for start & goal */ 
int random_seed; 
int density[16] [16]; 
struct node _ rec *node[66] [66]; 
float *ticks; 
{ 
struct node ree *“eurrent, “next, ceart, gear, 
int: k= 1; 


long sec, usec; 

Static struct timeval *tvp; 
Static struct timezone *tzp; 
Static Int fi rse, =. 


#ifdef IRIS 
union chrom_union dummy cu; 
dummy _cu.alleles = 0; 
#endif 


if ( first 
tvp = (struct timeval *)malloc(sizeof (struct timeval)); 
tzp = (struct timezone *)malloc(sizeof (struct timezone) ); 
first = F; 


srandom(random_seed); /* seed the random generator */ 
make _array(density, node); 

start = find_node(node, sx, sy); 

Start =-73b dee moron 

start->d1St Fromustare. sme, 
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goal = find_node(node, gx, gy); 
goal->state = GOAL; 

current = start; 

current->state = CURRENT; 


#ifdef IRIS 
/* initialize the IRIS system */ 
initialize (HCTITLE) ; 

#endif 


gettimeofday(tvp, tzp); 
sec = tvp->tv_sec; 
msec — Cvp->tv_usec; 


while ( current != goal) { 
next = move adjacent( node, current, start, goal ); 


if ( next != NULL ) 
next->predecessor = current; 
else if ( current->predecessor != NULL ) 
next = current->predecessor; 
else . 
printf ("\nNO SOLUTION - hill climb search"); 


Mext=>dist from start = 
CuErent—-dist stromestart teupdate dist ( current; next ); 

current->state = VISITED; 

current = next; 

current->state = CURRENT; 
#ifdef IRIS 

draw_ terrain (node, start,goal, current, next->dist_from_start, 

dummy _cu.factor); 
#endif 
} 


gettimeofday(tvp, tzp); 
*ticks = (float) (tvp->tv_sec - sec) + (tvp->tv_usec - usec) /1000000.0; 


#ifdef IRIS 
sleep ( SLEEPTIME ); 
#endif 


goal->state = GOAL; 
Sstart->state = START; 


Pecurm(sgoal—>dist from start ); 
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/* KKKKKKKKKKKKKKKKK move adjacent KKK KKK KKK KKK KK KKK KEKKKKKKKKEKKKKEKK * / 


struct node_rec *move_adjacent( node, c, start, goal ) 
struct node_rec *node[66] [66]; 
struct node _rec *c; /* Curreninw s) 
SEDUCE nOdenre> start, 
struct node rec *goal; 
{ 
Struce node Tec best, 
int xi, yi, base xi, base yi top x1, tepay= 


best = NULL; 

base x1 = c->xi 
base _ yi C=> V2 
Opes c=->xi 
top _yi C->yi 


1: c->xi - 
1: c=>yl - 
64 : c->xi + 
64 © °Cc=>yi + 


1 
i. 
64 
64 


e ‘“e eo 


ll 
ll 
il 
JV ev ov 
a ee eo 


‘we 


for (xi=base_xi;xi<=top_ xi; xitt) 
for (yi=base yi;yi<=top yi;yitt) 
Switch ( node[xi] [yi]->state ) 
{ 
case UNTOUCHED: 
case GOAL: 
node [xi] [yi] ->state = FRONTIER; 
node [xi] [yi]->dist_from_goal = update dist( node[xi][yi], goal ); 
best = find _best( best, node[xi] [yi] ); 
break; 
case FRONTIER: 
best = find _best( best, node[xi] [yi] ); 
break; 
} 


return( best );} 


/* KEKKKKKEKKKKKKKEKK fing best KKK KKK KKK KKK KK KKK KKK KKK KEKEKEKKKKKKEKKKKK * / 


/* Assigns the input node to best if appropriate */ 


Struct node_rec *find_best( best, n ) 
Struct Node recs eee 
SCEUCE Mode sre emai, 
{ 
if ( best == NULL } 
bese sn, 
else if ( n->dist_from_goal < best->dist_from_goal ) 
best saan, 
return( best ); 
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rt_astar.c 


/* 
File: 
Programmer: 
Environment: 
Language: 
Date: 
Revised: 
Comments: 


«/ 


Bemactar.c 
Gao] parker 
any 
Gc 
20 feb 92 

2 apr 92 


RTA-star search - Finds best adjacent node visited or frontier 


#include "ga_search.h" 
#include <sys/time.h> 


int heap_size; 


/* KKKKKKEKKKEKKKK KKK Yrta star KKK KKK KKK KKK KKK KKK KKK KKK KKK KKKKKKKKKK x / 


float rta_star( sx, sy, gx, gy, random_seed, density, node, ticks ) 
iit SsX,3Y,9X,gy; /* position in density array for start & goal */ 
int random_seed; 

int density[16] (16]; 

Seuuet modewrec *node[66] (66); 


itioat *ticks; 
{ 


weGuectenoge rec *Current, *Start, *goal; 
weoueceenode rec “best two[2]|; 


int k = l; 


long sec, 


usec; 


Static struct timeval *tvp; 
Seadtic Struct timezone *tzp; 
Seatic int first = T; 


#ifdef IRIS 


union chrom_union dummy_cu; 


dummy _cu.alleles = 0; 
#endif 
Pea first, )  { 
tvp = (struct timeval *)malloc(sizeof (struct timeval)); 
tzp = (struct timezone *)malloc (sizeof (struct timezone) ); 
first = F; 


heap_size 


= Q; 


srandom(random_seed); /* seed the random generator */ 
make array(density, node); 
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start = find _node(node, sx, sy); 
start->state = START; 
Start=-dise  tsemes tan cero 
goal = find_node(node, gx, gy); 
goal->state = GOAL; 

current = start; 

current->state = CURRENT; 


#ifdef IRIS 
/* initialize the IRIS system */ 
initialize (RTASTIZLE)- 

#endif 


gettimeofday(tvp, tzp); 
SGC = CvVp=-tv sce, 
usec = tCvp--tv usec, 


while ( current != goal) { 
best_two[0] = NULL; 
best_two[(1l] = NULL; 
update_rtastar_ adjacent( node, current, best two, start, goal ); 


if( best_two[0] == NULL ) break; 
if( best_two[1] == NULL ) 

current->dist from_goal = BIG NUMBER; 
efcoe 


current->dist_from_goal = best_two[1]->subtotal; 
current->state = VISITED; 
best_two[0]->dist_from_start = 
current->dist_from_start + update dist (current,best_two([0]) 
current = best two([0]; 
CUurrents=7state = CURRENT, 


#ifdef IRIS 
draw_terrain (node, start, goal, current,best_two[0]->dist_from_start, 


dummy cu.factor); 
#endif 


} 


gettimeofday(tvp, tzp); 
*ticks = (float) (tvp->tv_sec - sec) + (tvp->tv_usec - usec) /1000000.0; 


#ifdef IRIS 
sleep( SLEEPTIME ); 
#endif 


goal->state = GOAL; 


Start->state = START; 
return( goal—>dist from Stari, 
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[xk KkR RR KK KK KKK KK & update rtastar adjacent Rk KKK KK KK KK KK Ke / 


update _rtastar_adjacent( node, c, best two, start, goal ) 


Secuct 
Se ruct 
Seruct 
Seruct 
re uct 


{ 


node _rec *node(66] (66]; 

node_rec *c; /* current */ 
node _ rec *best_two[2]; 

foaerrec ~Start; 

node rec *goal; 


mitaxil, yi, base xi, base_yi, top xi, top yi; 


base x1 = c->xi == 1? #1: c=->xi - 1; 
base yi = c->yi == 172? 1: c->yi - 1; 
top x1 = c->xi == 64 ? 64 : c->xi + 1; 
top yi = c->yi == 64 ? 64 : cH->yi + 1; 


for (xi=base_ xi;xi<=top xi;xitt) 
for (yi=base yi;yi<=top yi; yitt) 


Switch ( node[xi] [yi]->state ) 


{ 
case UNTOUCHED: 
case GOAL: 
node [xi] [yi]->state = FRONTIER; 
node [xij [yi]->predecessor = c; 
node [xij] [yiJ->dist_from_goal = update _dist( node[xiJ[yi], goal ); 
node [xij [yi]->subtotal = 
node [x1] [yi]->dist_from_goal + update dist( node[xi] [yi], c ); 
insert ( best_two, node[xi] [yi] ); 
break; 
case FRONTIER: 
case VISITED: 
node [xi] (yi]->subtotal = 
Model (vary --dast trom goal + update dist ( node[xij[yi], ¢ ); 
insert( best two, node[xi] [yi] ); 
break; 
case CURRENT: 
break; 
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/* KKK KKKKEKKEKKKEK insert KKK KKK KKK KKK KKK KKK KEK KKK KKK KKK KKK KKKKKKKK * / 


/* Assigns the input node to best or second best as appropriate */ 


insert ( best_two, n ) 
struct node _rec *best_two[2]; 
struct node rec *n; 
{ 
if ( best_two[0] == NULL } 
best_two[0] = n; 
else if ( n->subtotal < best_two[0]->subtotal ) { 
best two[l] = best_two[0]; 
best two[0] =n; 
} 
else if ( best_two[1] == NULL ) 
best _two[l] = n; 
else if ( n->subtotal < best _two[1]->subtotal ) 
best two[l] = n; 
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bfsearch.c 


/* 
File: bfsearch.c 
Programmer: g.b. parker 
Environment: any 


Language: C 

Date: S july 92 

Revised: 

Comments: Best First Search - modified standard bfs for use with real- 
time search 
e/ 


#include "ga_search.h" 
#include <sys/time.h> 


/* Kwek KKK Ke KKK KKK KEKE KK bfsearch KEK KEKE KKK KKK KEK KKK KKK KKK KKK KKK KKK KKK KEK KKK x / 


float bfsearch( sx, sy, gx, gy, ind_chrom_factor, random_seed, density, node, 
ELcks ) 
Ant SxX,SY,9X,Qy; /* position in density array for start & goal */ 
Seebet tactor Struct ind chrom factor; 
int random_seed; 
int density([(16] [16]; 
Struct node rec *node[66] [66]; 
moat *ticks; 
{ 
SEruct node rec *current, *previous, *start, *goal; 
struct node_rec *frontier heap[4096]; 


long sec, usec; 

Static struct timeval *tvp; 
Static struct timezone *tzp; 
Beatie int first = T; 


float dist traveled = 0.0; 
eerie 7 


#ifdef IRIS 
union chrom_union dummy cu; 
dummy cu.alleles = 0; 


#endif 
met furst ) { 
tvp = (struct timeval *)malloc(sizeof (struct timeval)); 
tzp = (struct timezone *)malloc(sizeof (struct timezone) ); 
12 alesse, ES )a 
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Srandom(random_seed); /* seed the random generator */ 
heap_size = 0; 
make_array(density, node); 


Start = find_node(node, sx, sy); 
start->state = START; 
Start~>dist, fromistart = 0n0, 
goal = find_node (node, gx, gy); 
goal->state = GOAL; 


for (k=0;k<4096;k++) 
frontier heap[(k] = NULL; 
current = start; 
previous = start; 
#ifdef IRIS 
/* initialize the IRIS system */ 
initialize (BFTITLE) ; 
#endif 
gettimeofday(tvp, tzp); 
sec = tvp->tv_sec; 
usec = tvp->tv_usec; 


while ( current != goal) { 
bf update frontier list (node, current, previous, frontier heap, start, goal)% 
if( heap size == 0) { 
printf ("\nENDING SEARCH BEFORE GOAL - no more frontier"); 
break; 
} 
previous = current; 
current = bf pick_best_frontier(node,current,frontier heap,goal); 
dist_traveled = dist_traveled + current->dist_from_current; 
HLEGEE ERTS 
draw terrain (node, start,goal,current,dist_traveled, dummy _cu.factor); 
#endif 
} 


gettimeofday(tvp, tzp); 
*ticks = (float) (tvp->tv_sec - sec) + (tvp->tv_usec - usec) /1000000.0; 


#ifdef IRIS 
sleep( SLEEPTIME ); 
#endif 


goal->state = GOAL; 
start->state = START; 


return( dist _traveled ); 





[x KKK KK KKK KKKKKKKK HDF update frontier List *“** RRR KKK KKK KKKKKKKKKKEK x / 


Semupdeate frontier list( node, c, p, frontier heap, start, goal ) 
Beruct node rec *node([66) [66]; 
Beeuct node rec *c; /- Curccent 27 
Seruct node rec *p; /* previous */ 
meruct node rec *frontier heap[4096]; 
memuct node rec *start; 
Struct node_rec *goal; 
{ 
Paieesiyeeyt, base xi, base_yi, top_xi, top yi; 
float old_subtotal; 


Baseesti = C->xi - 1; 
base yl = c->yi - 1; 
Peemet — C->x1 + 1; 
Perey. = c->yi + 1; 


for (xi=base_xi;xi<=top_xi;xitt) 
for (yi=base_ yi;yi<=top yi;yitt) 
if ( (node[xi] [yi]->state == UNTOUCHED) || (node[xi] [yi] ->state == GOAL) 


node [xi] [yi]->state = FRONTIER; 

node [xi] [yi] ->predecessor = c; 

node [xi] [y1]->dist_from_goal = update_dist( node[xi] [yi], goal ); 
node [xi] [yi] ->subtotal = node[xi] [yi] ->dist_from_goal; 
insert_heap( frontier_heap, node[xi] [yi] ); 


/* RR K KK KK KEKEKKEKKEKKEK Ba. pick best €£rontier KKK KKK KKK KKK KKK KK KK KEKKKKEK x / 


/* Finds best frontier node, returns it */ 


Struct node rec *bf pick best frontier( node, current, frontier_heap, goal ) 
seruete node rec *node[66] [66); 
Semuersnoae rec *current; 
Struct node_rec *frontier heap[4096]; 
struct node _rec *goal; 
{ 
Seructemmode rec “best ptr, *q, *“gend, *qreset; 
Evoaeenede cost = BIG NUMBER; 
fioatenorm — 1.0; 9/* (current->dist_from_goal / 16.0); normalize factor */ 
Meek, Vi, ks 
float steps; 
int done = F; 
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best_ptr = frontier_heap[0]; 

q = Current, 

qend = current; 

QEreCsSee —seurrec ue, 
current —>back tudekustate = VIsTIEep, 
current->dist from_current = 0.0; 


while ( !done && (q '!'= NULL) ) { 


steps = q-=>dist from current 7) 1.0, 
for (k=0;k<8;k++) { 
if (k==4) 
steps = q->dist_from_ current + SQRT2; 
Xi = gen _xi( k, q->xi ); 
yi = gen yi( kj iq=-yaa): 


if( ( (node[xi] [yi]->state == VISITED) || 
(node [xi] [yi] ->state == FRONTIER) | | 
(node [xi] [yi]->state == START) ) && 
node [xi] [yi]->back_track_ state == UNTOUCHED ) { 
node [xi] [yi]->dist_from_current = steps; 


node [xi] [yi]->back_track_state = VISITED; 
node[xi] [yi]->qreset = qreset; 


qreset = node[xi] [yi]; 
if ( (node[xi][yi]->state == VISITED) || 
(node [xi] [yi]->state == START) ) { 


qend->qnext = node[xi] [yi]; 
qend = node[xi] [yi]; 
} 


else { /* node[xi] [yi]->state == FRONTIER */ 
if ( node[xi] [yi] == best _ptr ) 
done = T; 


} 


} 
q = q->qnext; 
} 
reset_back_ track_state( qreset ); 
best ptr->state = VISITED; 
delete heap( frontier heap, best _ptr ); 


EBCCUEN( DeE€sSERPEE a, 


psearch.c 


/* 
File: psearch.c 
Programmer: g.b. parker 
Environment: any 


Language: C 

Date: 9 july 92 

Revised: 

Comments: Persistence Search - uses distance to goal and distance to cur- 
rent 
to determine best frontier. 
B / 


#include “ga_search.h" 
#include <sys/time.h> 


/* KREKKKEKKKEKKKKKKKKK psearch KKK KKK KKK KKK KKK KK KK KKK KKK KKK KKK KKKKKKKKK * ff 


Bloat psearch({ sx, sy, 9X, gy, ind_chrom factor, random_seed, density, node, 
ticks ) 
init Sx,SV¥,gQx,9y; /* position in density array for start & goal */ 
Meedce factor struct ind chrom factor; 
int random_seed; 
int density([(16) (16); 
Beruct node rec *node[66] [06]; 
eloat *ticks; 
{ 
Struct node rec *current, “previous, *start, *goal; 
struct node rec *frontier heap[4096]; 


long sec, usec; 

static struct timeval *tvp; 
static struct timezone *tzp; 
Sicatac ant first = T; 

float dist_traveled = 0.0; 
ime K> 


me first ) { 


tvp = (struct timeval *)malloc(sizeof(struct timeval)); 
tzp = (struct timezone *)malloc(sizeof(struct timezone) ); 
fipstwF; 

} 

srandom(random_seed); /* seed the random generator */ 


heap size = 0; 
make_array(density, node); 

Start = find node(node, sx, sy); 
Seartsestate = START; 
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start=>disStatrom start o .o, 
goal = find_node(node, gx, gy); 
goal->state = GOAL; 


for (k=0;k<4096;k++) 
frontier heap[k] = NULL; 
current = start; 
previous = start; 
#ifdef IRIS 
/* initialize the IRIS system */ 
initialize (PERSTITLE) ; 
#endif 
gettimeofday(tvp, tzp); 
sec = tvp->tv_sec; 
usec = tvp tv usec, 


while ( current != goal) { 

if( adjacent (current, goal) ) { 
dist_traveled = dist_traveled + update _dist( current, goal ); 
current = goal; 

} 

else { 
p_update frontier list( node, current, previous, frontier_heap, start, 

Godly, Indechrommractoum ., 
if( heap size == 0) { 
printf("\nENDING SEARCH BEFORE GOAL %@d %d - no more frontier", 
current->xl, Current->yi); 
break; 
} 
previous = current; 
current = p_ pick best frontier( node, current, frontier heap, 
goal, ind chromsctaeror) 

dist_traveled = dist_traveled + current—>dist frem current; 

} 

#ifdef IRIS 
draw_terrain (node, start, goal, current,dist_traveled, ind _chrom_factor); 
#endif 
} 


gettimeofday(tvp, tzp); 

*ticks = (float) (tvp->tv_sec - sec) + (tvp->tv_usec - usec) /1000000.0; 
#ifdef IRIS 

sleep( SLEEPTIME ); 
#endif 

goal->state = GOAL; 

Start->state = START; 


return( dist_traveled ); 





/* KKK KKK KKKKKKK KKK p update frontier lee te KKK KKK KKKKKKKKKKKK KKK KKK x / 


mmupaace frontier list( node, c, p, frontier heap, start, goal, factor ) 
struct node rec *node[66]} [66]; 
struct node _rec *c; /  ecurrent => /; 
Struct node_ rec *p; /* previous */ 
Struct node_rec *frontier heap[4096]; 
memuce node rec *start; 
Struct node _ rec *goal; 
meee factor struct factor; 
{ 
Dieses, yi, base xi, base yi, top xi, top yi; 
eroac, Old subtotal; 


Base X1 “= ¢C=>x1 - 1; 
Basesyi = c->yli - 1; 
Eopexi = C=->x1 + 1; 
Pommyl = c=->yi + 1; 


for (xi=base_xi;x1i<=top xi;xitt) 
for (yi=base_yi;yi<=top yi;yitt) 
if ((node[xi] [yi] ->state == UNTOUCHED) || (node[xi] [yi]->state == GOAL)) { 
node [xi] [yi]->state = FRONTIER; 
node [xi] [yi] ->predecessor = ¢c; 
node [xi] [yi]->dist_from_goal = update dist( node[xi] [yi], goal ); 
node (xiJ [yi] ->subtotal = 
nece (xi) yt) --d1stmerom Goal = stactor. goal dist; 
insert _heap( frontier _heap, node[xi] [yi] ); 


/* KKK KK KKK KK KKK KKK p_ pick best frontier KKeKKKKKKKKKKKKKKKKKKKKK KKK K x | 


mamrmnds best frontier node, returns it */ 


struct node rec *p pick best frontier( node, current, frontier heap, goal, fac- 
tor ) 
Struct node rec *node[66] [66]; 
peruct modesrec *Current; 
struct node_rec *frontier heap[4096]; 
Peruectenode rec “goal; 
meruect factor Struct factor; 
{ 
PerNceenode “eCe=beSt ptr, *q, “qend, “qreset; 
Elode node Cost = BIG NUMBER; 
float norm = 1.0; /* (current->dist_from_goal / 16.0); normalize factor */ 
float lower bound = frontier _heap[0]->subtotal; 
float upper _bound = BIG NUMBER; 
Merri, Yi, ik; 
float steps; 
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DéSt Pr Ee —secUmrethe, 

q = current; 

gend = current; 

Greseee serene, 
CUrrent—>back track Stateu— Vio tik, 
current->dist_from_current = 0.90; 


while ( (lower_bound < upper bound) && (q != NULL) ) { 
steps = q->dist_from_current + 1.0; 
for (k=0;k<8;k++) { 
if (k==4) 
steps = q->dist_from_current + SORT2; 
xX1 = gen_xi( k, q->xi ); 
yi = gen_yi( k, q->yi ); 


if ( ( (node[xi] [yi]->state == VISITED) || 
(node[xi] [yi]->state == FRONTIER) || 
(node [xi] [yi]->state == START) ) && 
node [xi] [yi] ->back_track state == UNTOUCHED ) { 


node [xi] [yi]->dist_from_current = steps; 
node [xi] [yi] ->back_track_state = VISITED; 
node [xi] [yi] ->qreset = qreset; 

qreset = node[xi] [yi]; 


if ( (node[xi] [yi] ->state == VISITED) || 
(node[xi] [yi]->state == START) ) { 

qend->qnext = node[x1i] [yi]; 

gend = node[xi] [yi]; 
} 
else { /* node[xi] [yi]->state == FRONTIER */ 

node _ cost = node[xi] [yi]->subtotal + 

node [xi] [yi]->dist_from current * factor. currenteara, 


if (node_cost < upper_bound) { 
upper bound = node cost; 
best ptr = node[xi] [yi]; 

} 


} 
q = q->qnext; 
lower bound = frontier heap[0]->subtotal + steps * factor.current_dist; 
} 
reset_back track_state( qreset ); 
best_ptr->state = VISITED; 
delete _heap( frontier _heap, best_ptr ); 


met urn ( bestapres). 





tsearch.c 


/* 


Baile: tsearch.c 
Programmer: g.b. parker 
Environment: any 


Language: S 

Date: 9 july 92 

Revised: 

Comments: This is a multi-heuristic search that takes in the bias 


factors in the form of an eight digit hexadecimal number. 


/ 


#include “ga_search.h" 
#include <sys/time.h> 


int heap size; 


/* kkk kkk kkk kk kkk kkk search 


KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK x / 


float search( sx, sy, gx, gy, ind_chrom_factor, random_seed, density, node, 
Backs ) 


mit Sx,SY,9X,9Y; 


mecuat factor struct ind chrom factor; 
int random_seed; 

mat density[i6] [16]; 

SBeruct node rec *node([66]) [66]; 

mroat “ticks; 


{ 


struct node _ rec *current, *previous, *start, 
Beruet mede rec *frontier heap[4096); 


long sec, usec; 

Static Struct timeval *tvp; 
Static struct timezone *tzp; 
Stdcic int first = T; 


float dist_traveled = 0.0; 
int k; 


mf ( first) { 


/* position in density array for start & goal */ 


Sgoal; 


tvp = (struct timeval *)malloc(sizeof (struct timeval)); 
tzp = (struct timezone *)malloc(sizeof(struct timezone) ); 


first. = F: 
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srandom(random_seed); /* seed the random generator */ 
heap size = 0; 
make array(density, node); 


Start = find_node(node, sx, sy); 
Start->state = START; 
Sstart->dist from_start = 0.0; 
goal = find_node(node, gx, gy); 
goal->state = GOAL; 


for (k=0;k<4096;k++) 
frontier heap[k] = NULL; 

CULITENt,— stare, 

previous = start; 


#ifdef IRIS 
/* initialize the IRIS system */ 
initialize (GATITLE) ; 

#endif 


gettimeofday(tvp, tzp); 
Se€C = UvVp=7Ctv (sec, 
uSseC =" CVYp=-cy Usec, 


while ( current '!= goal) { 

if( adjacent(current, goal) ) { 
dist_traveled = dist_traveled + update_dist( current, goal ); 
current.= goal; 

} 

else { 
update frontier list( node, current, previous, frontier heap, start, 

goal, ind_chrom_factor ); 


if( heap size == 0) [{ 
printf ("\nENDING SEARCH BEFORE GOAL - no more frontier"); 
break; 

} 

previous = current; 


current = pick _best_frontier( node, current, frontier heap, goal, 
Ind chrom FAaCeoEs), 
dist traveled = dist traveled + current->dist_from_current; 
} 
#ifdef IRIS 
draw_terrain (node, start, goal, current,dist_traveled,ind_chrom_factor) ; 
fendif 
} /* end while loop */ 


gettimeofday(tvp, tzp); 
*ticks = (float) (tvp->tv_sec - sec) + (tvp->tv_usec - usec) /1000000.0; 
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#ifdef IRIS 
sleep( SLEEPTIME ); 
#endif 


goal->state = GOAL; 
Start->state = START; 


#ifdef SUN 
/* Print to standard output */ 
/* not normally used, but optional for sun */ 


/* 
meintt(™\n") ; 
print_node (node) ; 
Beinet ("\n") ; 
Perott({"\nDIST = %£", dist traveled); 
:/ 
#endif 


return( dist traveled ); 
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tfrontier.c 


/* 
File: tirontierge 
Programmer: g.b. parker 
Environment: any 


Language: & 

Date: 9 july 92 

Revised: 

Comments: Maintenance of frontier list 
| 


#include “ga_search.h" 
#include "math.h" 


/* KKK KK KKK KKKKKKKK update frontier list KKK KKK KK KK KKK KKK KKK KKEKKKKKK 


/* Looks two away from the current node to update stable search 
characteristics */ 


update EFronticr list (@ncde; cp, frontier heap, start, goal, factor } 
Struct node _ rec *node[66] [66]; 
struct node” rec “cr; /* current */ 
struct node_rec *p; /* previous */ 
Struct node_rec *frontier heap[4096]; 
SEXTUCGCENnOde” ree cteame, 
SEFUCE Mode Tecmo gear, 
Struct £aCtorStruce factor, 
{ 
int Xi, yi, base xi; base yi topixt,wetope, =, 
float old subtotal; 


1: c->xi - 
aT 1 : c=>yl - 
64 64 : c->xi + 
64°27 64 3c ->yr 


Ih 


J) 


base xi = C->xi 
base yi = c->yi 
topes G=exi 
top yi s= "e=7ya 


Il 

I 

ll 
Oe) 


update_adjacent_obstacles( node, c ); 
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Dy 


for (xi=base_xi;xi<=top_xi;xitt) { 
for (yi=base yi;yi<=top_yi;yit++) { 


if ( ((node[xi] [yi]->state == UNTOUCHED) | | 
(node [xi] [yi] ->state == GOAL) ) && adjacent (node[xi] [yi],c) ){ 

node [xi] [y1]->state = FRONTIER; 

node [xi] [yi] ->predecessor = c; 

node [xi] [yi]->dist_from_goal = update dist( node[xi] [yi], goal ); 

node[xi] [yi]->dist_from_start = update _dist( node[xi] [yi], start ); 

node [xi] [yi]->subtotal = 

compute_subtotal( node[xi] [yi], factor, 
calc_momentum(node[xi] [yi],c,p), 
update crowd sides (node,node{x1i] [yi]), 

update crowd diag(node,node[xi][yi]) ); 

insert_heap( frontier heap, node[xi] [yi] ); 


else if ( node[xi] [yi] ->state == FRONTIER ) { 
Sldesubcotal = node[xi)] [yi]->subtotal; 
node [xi] [yi] ->subtotal = 
compute subtotal( node[xi] [yi], factor, 

calc _momentum(node[xi] [yi],c,p), 
update crowd sides (node, node[xi] [yi]), 
update crowd _diag(node, node[xi] [yi]) ); 

if ( !tequalf(old_subtotal,node[xi] [yi]->subtotal) ) { 

move heap( frontier heap, node[xi] [yi]->frontier_index )j; 


} 
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/* KKK KKKKKKK KKK KK KK update _adjacent_obstacles RK KK KKK KKK KK KKK K x / 
/* Records for future use which adjacent nodes are obstacles */ 


update_adjacent_obstacles( node, c ) 
struct node_rec *node[66) [66]; 
struct node_rec *c; /* current */ 
{ 
int xi, yi, base_xi, base_yi, top xi, top_yi; 


base xi > ¢—>xi ae, 
base. yi.) Cl>yig Se 
top. Xi ) =fe—> xi, 
top yi  ~=oe--yae 1; 


for (xi=base x27 x1<=top xi; sci +) 
for (yi=sbase_yi-syi<=top_yi;yitt) 
if ( node[xi) [yi]->state == OBSTACLE ) 
node [xi] [yi] ->back_track_state = OBSTACLE; 


/* KEKE KKK KK KK KK KK pick best frontier KKK KKK KKK KKK KKK KK KKK KKK KKKKKK x / 


/* finds best frontier node, returns it */ 


Struct node _ rec *pick_best_frontier(node,current, frontier heap,goal, factor) 

struct node_rec *node[66] [66]; 

StLUCE Nedemrec. =cucEcnc, 

Struct node_rec *frontier_ heap[4096]; 

Struct nodé rec -=goal, 

SErUGCE Lactorvsenuce factor, 

{ 
Struct nodesrec “best ptr, “Gq, eadend.s .qreser, 
float node_cost = BIG NUMBER; 
float norm = 1.0; /* (current=>dist from goal / 16°70); normalize factor */) 
float lower bound = frontier _heap[0)->subtotal; | 
float upper bound = BIG NUMBER; 
int xi, yoy 
float steps; 





best [ptry =sermerent, 

q = current; 

qend = current; 

qreset = current; 
CUETeENT =7back track Sstaeeg= =) Foti eD, 
current->dist_from_current = 0.0; 
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while ( (lower bound < upper _bound) && (q != NULL) ) { 


Steps = q->dist_from_current + 1.0; 
for (k=0;k<8;k++) { 
if (k==4) 
steps = q->dist_from_current + SQRT2; 
Eemmgen xi( kK, G->xi ); 


yi = gen_yi( k, q->yi ); 
if ( ( (node[xi] (yi]->state == VISITED) || 
(node [xi] [yi] ->state == FRONTIER) || 
(node [xi] [yi])->state == START) ) && 
node [xi] [yi] ->back_track_state == UNTOUCHED ) { 
node[xi] [yi]->dist_from_current = steps; | 


node [xi] [yi] ->back_track_state = VISITED; 
node [xi] [yi]->qreset = qreset; 
qreset = node[xi] [yi]; 


if ( (node[xi] [yi]->state == VISITED) || 
(node [xi] [yi]->state == START) ) { 
qend->qnext = node[xi) [yi]; 
qend = node[xi] [yi]; 


} 
else { /* node[xi]) [yi]->state == FRONTIER */ 


node cost = node[x1i] [yi]->subtotal + 
node [xi] [yi]->dist from_current * factor.current dist + 
calc_move_away (node[xi] [yi],current,goal) * factor.move away; 
if (node_cost < upper bound) { 
HpeeEseound — nede cost; 
best ptr = node[x1] [yi]; 
} /* end if */ 
} /* end else */ 
} /* end if */ 
} /* end for loop */ 
q = q->qnext; 
lower bound = frontier heap[(0]->subtotal + steps * factor.current_dist; 
} /* end while loop */ 
reset Dack track state( qreset ); 
best _ptr->state = VISITED; 
delete heap( frontier_heap, best ptr ); 


Beturn( DESt per); 
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/* KKK KKK KKKKAKKEK reset back track state KKK KKKKK KKK KKK KKK KK KKK KEK x / 


/* Resets node record fields used to perform the backtrack search */ 


reset back  trackescare( qrescers) 
SCrUCE NcdeRLrec “Gresser, 
{ 


struct node rec *temp; 


while (qreset != NULL) { 
temp = qreset->qreset; 
Greset=>back track State = UNTOUC EE, 
qreset->qreset = NULL; 
qreset->qnext = NULL; 
qreset = temp; 





Via Kaa KK KKK KKK KK KKK update list Ket KK aK KKK KKK KKK KK KKK KKK KKK KKK KKK KK x / 


/* Euclidean distance between input nodes */ 


float updatesdise( ni; en2)) 
struct node rec *nl; 
struct node _rec *n2; 


{ 


nl->x - n2->x; 
nls) eZ 


floaty x 
float y 


return( sqrt ( x*% + y*yv oe 


/* KKK KK KK KKK KKK KKK roles Tote He KK KK a a KKK KKK KKK KKK KKK KKK KKKKKKKKKKKK * / 


/* Absolute difference between two integers */ 


diff Vint Cua} 
int fal, b; 
{ 


int c =a -b; 


if (e3< 70») 
return (s-ce) 
else 
return( c ); 
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/* KKK KKK KKK KKK KKK adjacent KKK KKK HK KKK KKK KKKKAKK KK KKK KKK KKK KK KK x / 


/* Returns T if the input nodes are adjacent */ 


adjacent (nl,n2) 
meeuct node rec *nl; 
meruct node rec *nZ; 


{ 
Boe (diet ant (nl—>xi,nZ2—->xi) < 2) && (diff int (nl—->yi,nZ2<—>yi) < 2) ); 


[& kk RK Update Crowd Sides RRIF III III III III TICK II ae / 


/* Counts the known adjacent horizontal/vertical obstacles to the frontier 
node */ 


update crowd sides(node, f) 
Struct node rec *node[66] [66]; 
mepuct node rec *f; /* frontier */ 


{ 


int s_count = 0; 

oe * / 

if ( node[f->x1i] [f->yi + 1])->back_track_ state == OBSTACLE ) 
eecount = s count + 1; 

oe. */ 

if ( node[f->xi + 1] ([£f->yi]->back_track state == OBSTACLE ) 
Smeounte= Ss count + 1; 

ie | 

ai4 node [f->x1] [f->yi - 1]->back track state == OBSTACLE ) 
ee-coune — S count + 1; 

aN * / 

iiaeenode|i—->xil — 1) (f->y1)—>back track state == OBSTACLE ) 


Scount = 3 count + 1; 
Beeurhn( Ss count ); 


[I IIIT update crowd diag %% % 5 III I I I I IOI IOI II Ie / 


/* Counts the known adjacent diagonal obstacles to the frontier node */ 


update_crowd_diag (node, f) 

struct node_rec *node[66] [66]; 
struct mode tec *f-an/. = © onete maa, 
{ 


int d_count = 0; 


[nN ae 

if ( node[f->xi + 1] [f->y1 + 1]->back_track_state == OBSTACLE ) 
dVecount = 7d Counts, 

[aS 7 

if ( node[f->xi + 1] [f->yi - 1]->back_track_state == OBSTACLE ) 
d count = dicount. 7 7, 

f= SWes/ 

if( node(f->xi - 1) [£->y1 = 1]=>back track ‘statce—— @5of-cuse, 
d counts = "GC coumt aa, 

{= NW, 


if( node[f->xi - 1] [(f->yi + 1]->back_track_state == OBSTACLE ) 
Ga count = dlcouneG stoi, 
return( d count ); 


/* KKEKKKKK KKK KKK KK KK calc move away KKEKKK KKK KKK KKK KKK KKK KE KKK KKK KK KK KK x / 


/* Determines if a move to the frontier would be moving away from the 
goal. Each axis move away counts as two. */ 


calc _ move _away( f, c, g) 
SCEUCE Boge mHCc ma, 
Struct node rec *c; 
S€rUCESNOde IreCuaG, 

{ 


8 


Int malcoune = 0; 

1f( diff int (f->xi,g->x1) 9 > ditieint (co -<' ge) 
ma_count = 2; 

if ( diff_int (f->yi,g->yi) > diff_int(c->yi,g->yi) ) 


ma_count = ma_count + 2; 


FECurEn | maTcount. 7 
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[xk kKkKKK KKK KKK ERK KKKK CalC MOMENCUM —* % RK III III OKO III KR a / 


/* Returns 0 if no change in direction, 1 if 45 degree change, two if 90 
degree change, and three if 135 degree change or node not adjacent */ 


calc_momentum( f, c, Pp ) 
Breruct node rec *f; 
menue: node rec *c; 
meruct node rec *p; 
{ 
if( adjacent(f,c) && adjacent(c,p) ) 
Peon oer nth x eee — oxi, C->xi —ef—->xi) + diff int(p->yi <- ¢c->yi,c- 
Byt = f->yi) ); 
else 
return( 3 ); 


/* REKKKKKKKKKKKRKKRKK compute subtotal KKKKKKKKKKKKRK KKK KKK KKK KKK KK KKK x / 


/* Computes the frontier nodes subtotal value dependent on the stable heuris- 
SrCcs e/ 


PlIademcompute: SuDbtotal( n, factor, m, cs, cd ) 


struct node rec *n; /* the node */ 

Seralce Lactor struct factor; 

Ent m; /* momentum */ 

mt Cs; /* crowding sides */ 

mnt cd; /* crowding diagonals */ 


{ 
Pewee N—--Gist from start * fLactor.start dist + 
irr me romugad len Lactor.goal Gist + 
CGsu* factor.crowa sides + 
ede~ factor. crowd diag + 
m * factor.momentum ); 
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theap.c 


/* 
File: theap.c 
Programmer: g.b. parker 
Environment: any 


Language: & 

Date: 9 july 92 

Revised: 

Comments: Frontier heap functions 
7 


#include "“ga_search.h" 


int heap_size; 


/* KKK KKK KKK KKK KK KK insert heap tem KKK KK KKK KKK KKK KKK KKK KKK EEK x / 


/* Inserts a node into the frontier heap */ 


insert_heap( fh, n ) 
struct node_rec *fh[(4096]; /* frontier heap */ 
struct node_rec *n; 7 /* node to insert */ 
{ 

n->state = FRONTIER; 

n->frontier_ index = heap size; 

fh{heap size] = n; 

heap _ size = heap_size + 1; 

move_heap( fh, heap_size-1l ); 


/* wR KKK KKK KKK eae vtere heap wok am mK am aK KK KKK KKK KKK KKK KKK KKK x / 


/* Deletes a node from the frontier heap */ 


delete heap( fh, n ) 
struct node_rec *fh[(4096]; /* frontier_heap */ 


struct nodegrec *n; /* node to delete */ 
{ 

heap _ size = heap size - 1; 

fh(n->frontier index] = fh{heap size]; 


fh(n->frontier_index]->frontier index = n->frontier_index; 
n->state = VISITED; 
fh(heap_ size] = NULL; 
if ( n->frontier_index != heap size ) 
move_heap( fh, n->frontier_index ); 
n->frontier_index = NA; 
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/* Kaw KKK KKK KKK KK move heap KA K KK KKK KKK KEK KKK KKK KKK KK KKK KKK KKK KKK x / 
/* Moves a node in the frontier heap if required; dependent on the value 
of the nodes subtotal field. */ 


move _heap( fh, i ) 
struct node rec *fh(4096]; /* frontier_heap */ 
mit i; /* index of node to possibly move */ 
{ 
int parent (A=) 2; 
Boe Cha ld G(2ritieee= heap Size)? 1: 2*i+1; 
imcechiitd2 a (2*i+2) >= heap size) ? i : 2*i+2; 
if ( (child2 != 1) && (fh[child2]->subtotal < fh[child]->subtotal) ) 
ehild = child2; 
while (fh[i]->subtotal < fh[parent]->subtotal) { 
swap( fh, i, parent ); 
i = parent; 
parent = (i - 1) / 2; 


} 
while (fh[iJ->subtotal > fh[child]->subtotal) { 


swap( fh, 1, child ); 


= child; 

ehaliadm = ((2*it1) >= heap size) ? i : 2*itl; 

etmeeta2 — {((2*1+2) >= heap size) ? 1 : 2*1i+2; 

femeechalds §!= a) && (fh[child2)->subtotal < fh[{child]->subtotal) ) 


child = child2; 


/* KKK KK KKK KKK KAKA KK swap KEKE KKK KKK KKK KKK KAKA KKK KKK KKK KKK KK KKK KKK KKK x / 


/* Swaps nodes in the frontier heap */ 


eeapt fh, il, i2 ) 
Beruct node rec *fh[(4096)% /* frontier heap */ 
mnt il, i2; /* indexes of nodes to swap */ 
{ 
Struct node_rec *temp ptr = fh[il]; 


I 
- 
No 


fh[il)->frontier index = 
eu(izZ}=>frontier index = il; 
eal) = thiliz): 

Cos 2) >= temp per; 
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evolve.c 


/* 
File: evolve.c 
Programmer: g.b. parker 
Environment: any 


Language: c 

Date: 9 july 92 

Revised: 

Comments: Performs selection, crossover, and mutation on input population. 
i! 


#include “ga_search.h”" 





/* kkkkkkkkkkkkkkkkk ayolve KK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK x / 


/* takes in a population, with fitness information which is used to produce 
the next population */ 


evolve( individual, rs ) 
Struct individual est ructa indivtcualio2 ie, 
ine. 6S; /* random number */ 
{ 
static struct individual structs *tempeina oa 
Static int Elest = 17 <9 tae ere 
Init. 
int, top = 0; 
i gNe, fevers, Cjeiole, i= Hus 
int mate _heap[31]; 
int even = T; 


PE ORG est es) ea 
create population( temp_ind ); 
Elest. = iE, | 

} 

temp_ind[0]->fit_sum = 

temp_ind[(1]->fit_sum 


il 
oe) 
Ooo 


create mate _heap( mate_heap, individual, rs ); 
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For (k=0;k<3Z;k++) { 
top = top + individual [k])->fitness; 
while ( (mate_heap[1] <= top) && (mate_heap[0] > 0) ) { 
temp ind(next_spot]->chrom.alleles = individual[k]->chrom.alleles; 
temp_ind[next_spot]->previous_index = k; 
if ( even ) { 
next spot = next_spot + 2; 
even = next_spot < 30 ? T : F; 
} 
else 
HeExt Speme— get odd); 
pop_mate_heap(mate_ heap) ; 
} 
pee 2ndavyidual([k]=>fit sum > temp ind[0])=>fit_ sum ) { 
set_equal( temp_ind[0], temp_ind[1], temp_ind[(0]->previous_ index ); 
set_equal( individual[k], temp _ind[0], k ); 
} 
else if ( (individual[k]->fit_sum > temp_ind[1]->fit_sum) && 
(individual[k]->chrom.alleles != temp_ind[0)->chrom.alleles) 
set_equal( individual[{k], temp _ind[1], k ); 
} 
crossover( individual, temp_ind, rs ); 
/* mutate done in crossover */ 


/* KKEKKKKKEKKKKKKKKKK create mate heap KKK KKK KKK KKK KK KKKKKKKKKKKKKEK KKK * / 


/* Creates a heap of integers which will be used to stochastically choose 
individuals for reproduction */ 


Create mate _heap( mh, ind, rs ) 


ent omna(S.1 ) ; 
BeGuctmrndividual “struct *ind[32]; 
nt 6S; /* random_seed */ 
{ 
int k; 


Pieweecal £it — OF 


srandom(rs); 
mh[(0]) = 0; 


for (k=0;k<32;k++) 
rota | ewe —storal fit + ind[{k])->fitness; 


fou (k= 2 ks) ey) 
insert_mate_heap( mh, (random() % total_fit) ); 
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) 


[RRR KK KKK KKK KKEKKKKKKEKKKEKEK Crossover FUNCCIONS KKK KKKKKKKKKKKKKKKKK KK KK KK / 


/* KEKE KKKKKKKKKKKK allele crossover kKaekkKkKK KKK KKK KK KKK KK KKK KKK KK KKK x / 


/* Performs crossover of the chromosome at a random allele position */ 


allele crossover( ind, temp_ind, k, cross allele ) 

struct individual struct “indieci, 

struct individualgstruct “Cempeminalo2y, 

int k, cross_allele; 

{ 
ind(k]->chrom.alleles = temp_ind[k]->chrom.alleles; 
ind[k+1]->chrom.alleles = temp_ind[k+1]->chrom.alleles; 


switch ( cross allele ) 
{ 
case 0: 
ind[k]->chrom.factor.start dist = temp_ind[k+1]/—>chrom. factor staresa umm 
ind [(k+1]->chrom.factor.start_dist = temp_ind[k]}=>chrom. factor. stared 
Caceeu. : 
ind[k]->chrom.factor.goal_ dist = temp_ind[k+1]—>chrom. factor. gedaleamee, 
ind[k+1]->chrom.factor.goal dist = temp_ind[k]—>chrom. factor Goaimar-c, 
Case Z- 
ind(k]-~>chrom.factor.current dist = 
temp ind[k+1]->chrom.factor.current_ dist; 
ind[(k+1]->chrom.factor.current dist = 
temp_ind[k]->chrom.factor.current_dist; 





case 3: 
ind[(k]->chrom.factor.crowd sides = 
temp ind[k+1]->chrom.factor.crowd_sides; 
ind[k+1]->chrom.factor.crowd_sides = 
temp_ind[k]->chrom.factor.crowd_sides; 
case 4: 
ind[k]->chrom.factor.crowd_diag = temp _ind[k+1]->chrom. factor.crowd_diag; 
ind(k+1]->chrom.factor.crowd_diag = temp_ind[k]->chrom.factor.crowd_diag; 
case 5: 
ind[(k]->chrom.factor.move away = temp _ind[k+1]->chrom.factor.move_away; 
ind(k+1]->chrom.factor.move_ away = temp_ind[k]->chrom.factor.move_away; 
case 6: 
ind[k]->chrom.factor.momentum = temp _ind[k+1]->chrom.factor.momentum; 
ind[(k+1]->chrom.factor.momentum = temp_ind[k]->chrom. factor.momentum; 
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/* KKK KKK KKK KK KKK bit crossover KKK KK KKK KKK KKK KKK KK KKK KKK KKKKKKKK x / 


/* Performs crossover of an allele at a random bit position */ 


bit crossover( ind, temp_ind, k, cross_allele ) 
meruct individualstruct *ind[32]; 

Beruct individual struct *temp_ind([32); 

Bit k, cross allele; 

{ 


Mmecross bit, inv cross bit; 


cross bit = get_mask( rand5({) ); 
inv cress bit = cross bit * MASK4; 


switch ( cross_allele ) 
{ 
ease. . 
maGiikl=echrom.factor.start dist = 
Vecempeind|k])—>chrom.factor.start dist & cross bit ) | 
(Beano eid (k+l) —-chrom-factor.startmdistws inv cross bat ); 
Poidiikt.)=—-Cnirom. factor .start dist = 
(stemp ind[k+1])->chrom.factor.start dist & cross_bit ) | 
(ecnpmincin)—-chiom.tactor.start@edist & inv cross bit ); 
break; 
Case 2: 
Back) =-Chrom.tfactor.goal dist = 
( temp_ind[({k]->chrom.factor.goal_ dist & cross bit ) | 
(eenomind(k+1)—-chrom.ftactor.goal dist & inv_cross bit ); 
ind({k+1]->chrom.factor.goal dist = 
( temp_ind[(k+1]->chrom.factor.goal_dist & cross_bit ) | 
ce mMomincit) —-Chrom factor goal. dist & inv crossmbit wy; 
break; 
Case 3: 
immeiKl|=-chrom.factor.current dist = 
mEcmpmria lk )->chrom. factor current dist & cross bit ) | 
(eecipeindik+il=>cheom factor-current dist & inv _cross_bit ); 
na (Kt) }i->chrom.facter.current dist = 
( temp_ind({k+1]->chrom.factor.current_dist & cross_bit )}) | 
(eremumtna (icl—->Cchrom.factor.current dist & inv cross bit )- 
break; 
case 4: 
ind[(k]->chrom.factor.crowd sides = 
cen ind k)=>chrom.factor.crowd sides & cross bit ) | 
ectpmticiikt | —-Chrom.factor.crowd Sides & inv _cross_bit ); 
ind(k+1]->chrom.factor.crowd sides = 
( temp_ind[(k+1]->chrom.factor.crowd_sides & cross bit ) | 
( temp_ind[(k]->chrom.factor.crowd_ sides & inv_cross bit ); 
break; 
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case 5: 
ind[(k]->chrom.factor.crowd diag = 
( temp ind[k]=>chrom. factor. crowd didg Smerosssolt) | 
( temp_ind[k+1]->chrom.factor.crowd_diag & inv_cross_bit ); 
indiik+1]->chrom.facver. crowd diag = 
( temp_ind[k+1]->chrom.factor.crowd_diag & cross_bit ) | 
( temp_ind[k]->chrom.factor.crowd_diag & inv_cross_bit ); 
break; 
case 6: 
ind[k]->chrom.factor.move_away = 
( temp_ind[k]->chrom.factor.move_away & cross _bit ) | 
( temp _ind[k+1]->chrom.factor.move_away & inv_cross_ bit ); 
ind[(k+1]->chrom.factor.move away = 
( temp_ind[k+1]->chrom.factor.move away & cross bit ) | 
( temp_ind[k]->chrom.factor.move_away & inv_cross_bit ); 
break; 
case 7: 
ind[(k]->chrom.factor.momentum = 
( temp_ind[k]->chrom.factor.momentum & cross bit ) | 
( temp ind[k+1]->chrom.factor.momentum & inv_cross_ bit ); 
ind[k+1]->chrom.factor.momentum = 
( temp_ind[k+1]->chrom.factor.momentum & cross bit ) | 
( temp ind[(k]->chrom.factor.momentum & inv_cross_bit ); 


/* waa KEKKKKKKKKKKKK Crossover KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KEKKKKKKK xy, 


/* The main function */ 


crossover( individual, temp_ind, rs ) 
struct individual struct *individuatie 2). 
struct individual struct “temp andi, 
int ws /* random_seed */ 
{ 

Int ky; 

int cross allele; 


srandom(rs) ; 
Cressmamrele = randsige, 


set _equal( temp_ind[0], individual[0], temp_ind[0]->previous_index ); 
set_equal( temp_ind[1], individual[1], temp_ind[1]->previous_index ); 


for (k=2;k<32;k=k+2) { 
allele crossover( individual, temp_ind, k, cross_allele ); 
bit crossover( individual, templind, Kk, crossed tele: 
mutate ( individual, k, rs ); 
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as KKK KKK KKK KKKKKKKEK get mask KKK KKK KKK KK KKK KKK KKK KK KKK KKK KKK KKK KKKKEK x / 


y* called by bit ‘crossover */ 


int get_mask( rn ) 
int rn; /* random number */ 
{ 
swatch ( rn ) 
{ 
case 0: 
return MASKO; 
case l: 
return MASKI1; 
case 2: 
return MASK2; 
case 3: 
return MASK3; 
case 4: 
return MASK4; 


/* KEK KKKKEKKKKKKKKKK mutate KKK KKK KKK KKK KKK KKK KKK KKK KKK KEK KKKKKKKKEKKKEK x / 


/* Runs through each bit of the chromosome determining if if will invert */ 


emutate ( ind, k, rs ) 

peoeect Individual struct *ind[(32]; 

mick, xs; 

{ 
“Unsigned int mut _factorl (Beet eeer £0: 
unsigned int mut_factor2 = Oxffffrffto; 
Pe gs 


for (g=0;g<28;g++) { 
mut_factorl = ( mut_factorl << 1) + one _i1f mutate(); 
Me eweadechOrcs— emut tacterz << 1 } + one if mutate(); 
} 
ind(k]->chrom.alleles = ind[(k]->chrom.alleles * mut_factorl1; 
ind(k+1]->chrom.alleles = ind[(k+1l1]->chrom.alleles %* mit faetcorZ; 
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/* KKK KKK KKK KEKE one ase mutate KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK x / 


/* Returns 1 if mutation is to take place at the present bit */ 


one_if mutate () 
{ 
if ( rand10000() < PROB BIT MUTATE ) 
return( 1 ); 
else 
return ( 0) > 


/* KKK KKK KKK KKK get odd KKK KKK KKK KKK KEK KKK KKK KKK KKK KK KKK KKK KK KK x / 


/* Determines placement of selected individual for reproduction. 
Distributes individuals to avoid mating of like chromosomes. */ 


/* Definitions only pertinent to this function */ 
#define LOW 0 
#define MED 1 
#define HIGH 2 


get_odd() 


{ 
static int next = LOW; 
Static int base 


Il 
bh 


Switch ( next ) 
{ 
case LOW: 
next = MED; 
base = base + 2; 
return( base ); 
case MED: 
next = HIGH; 
return( base + 10 ); 
case HIGH: 
next = LOW; 
if ( base == 11) { 
base = 1; 
return ( 31a): 
} 
else 
return( base + 20 ); 
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/* KKK KKK KKK KKKKKKEK set-equal KK KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKKKKK x / 


/* Sets one individual equal to another */ 


set_equal( from_ind, to_ind, k ) 

Bernet individual struct *from_ind, *to_ ind; 

mit k; 

{ 
Pomrid—-chrom.alleles = from _ind->chrom.alleles; 
to _ind->fitness = from_ind->fitness; 
Pomina->tit sum = from_ind->fit_ sum; 
to_ind->previous_index = k; 


Ill 


eheap.c 


/* 
File: eheap.c 
Programmer: g.b. parker 
Environment: any 


Language: C 

Date: 9 july 92 

Revised: 

Comments: Frontier heap functions 
if 


#include “ga_search.h" 


/* KKKKKKKKKKKKKKKEKK insert mate heap KKK KKK KK KKK KKK KKKKKKEKKKKKKK KKK x / 


insert_mate_heap( mh, num ) 
int, mn [sty /* matemneara,, 
int num; /* number to insert */ 
{ 
mn[0] = mh[0) + 1; 
mh{ mh[(0) ] = num; 
move mate heap( mh, mh[0] ); 
} 


/* KRKEKKKKKKKKK KK KKK K pop_mate heap KKK KK KKK KKK KKKEKKKKKKEKKKEKKKEKKKKKEKKK x / 


/* Removes top of mate heap */ 


pop_mate_heap( mh ) 
int mh[31]; /* mate_heap */ 
{ 

monf{l] = mh[{ mh[0] ]; 

mi{ mh[0] ] = 0; 

mh[(0] = mh[0] - 1; 

move _mate _heap( mh, 1 ); 
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[* kkk KKK KKK KKK IK Move mate heap  **— 5 — IIIT III III II IO ok / 


/* Readjusts heap after addition/removal of one of its members */ 


move mate heap( mh, i ) 


mnt mh[31); /* mate_heap */ 
mit i; /* index of num to possibly move */ 
{ 

miceparent =i =2= 1?1: i/ 2; 


mae child C(Z*1) >> mnl(0)]) 2? 1 = 2*i; 
tapeenitla2 = ((2%*1+1) > mh(0)) ? 1 : 2*1i+1; 


mmccecn1ld2 ‘= 1) &€& (mh[(child2) < mh(child]) ) 
eaela = child2; 


while (mh[i] < mh[parent]) { 
swap _num( mh, i, parent ); 
i = parent; 
Barent = 1== 171: i/ 2 


while (mh[{i] > mh[child])) { 
Swap _num( mh, 1, child ); 


ie Child; 

eariae = §((2*1) > mh(0)) ? 1 : 2*1; 

child2 = ((2*i+1) > mh(0)) ? i : 2*1i+1; 

if ( (child2 '!= i) && (mh[(child2] < mh[child]) ) 


ehaia = child2; 


yx KKK KKK KKK KKKKKKK Swap num KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK x / 


/* Swaps positions of two members of the mate heap */ 


Swap num( mh, il, i2 ) 
int mh[(31); /* mate_heap */ 
int il, i2; /* indexes of numbers to swap */ 
{ 
int temp_num; 


temp_num = mh[il]; 
mn{ilj) = mh[i2]; 
mh (i2] temp_num; 
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tmisc.c 


/* 
File: EMLSse .c 
Programmer: g.b. parker 
Environment: any 


Language: C 

Date: 6 apr 92 

Revised: 

Comments: Miscellaneous functions 
uid 


#include “ga_search.h" 


/* KKK KKKEKKEKKEKKKEKKEKEK gen XL RK KEK KKK KKK KKK KKK KKK KK KKK KEKE KEKE KEK KKK KKK KK x / 


/* Generates an integer value dependent on the input xi and k. 
Used with gen yi to generate all adjacent nodes to (xi,yi). */ 


int gen_xi( k, xi ) 
intr k; 
EnGe- xi; 
{ 
Switch (k) 
{ 
case 0: 
return (2x2) 
case l: 
return( xitl ); 
case 2: 
return<) x11)" 
case 3: 
return( xi-l ); 
case 4: 
patsnebaige((. pests oily jt 
case 5: 
return xi+1~—); 
case 6: 
return( xi-l ); 
case 7: 
return( xi-l ); 
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/* REAR KKAKKKAKKKKKEKK gen yi RRR KKK KKK KEKE KKK KKK KKK KKK KKK KKK aKa KKK x / 


/* Generates an integer value dependent on the input yi and k. 
Used with gen xi to generate all adjacent nodes to (xi,yi). */ 


int gen_yi( k, yi ) 
mot. k; 
mit yi; 
{ 
switch (k) 
{ 
case 0: 
eetwrn (hyitl ) ¢ 
case l: 
return ( yi ); 
case 2: 
Becueatyil—1) ; 
case 3: 
return( yi); 
case 4: 
retuin( yitl ); 
ease 5: 
return( yi-l ); 
case 6: 
Betwri( viol, ) >; 
case 7: 
reaurnir vat); 


/* KRaAaRKKKKaKKAKK KKK KKK equalf KKK KK KKK KKK KKK KKK KKK KKK KK KEK KKKKKEKKAKKKEEK EK x / 


/* Checks if two floats are equal (within 0.0001) */ 


int equalf( x, y ) 
bloat. x, Vv; 
{ 


momme( (x-v) < —0.0001) || ((x=-y) > 0.0001) ) 
return( F ) +; 
else 


meturn¢ TL );? 
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/* KKK KKK KKK KKK KK show least nodes KK KKK KKK KKK KKK KKK KKK KKK KKKKK « / 
/* Sets the state field to x for all nodes in the shortest path */ 
/* Not currently used, but available for graphics */ 


show_least_nodes( node, g ) 
struct node_rec *node[66] [66]; 
Struct node rec *g; /* goal */ 
{ 

struct node rec “best _ptz; 

float oes, 

int. s<isyviyk. 


while ((g->dist_from_start > 0.0) && (g->dist_from_start < 10000.0)) { 
best = BIG NUMBER; 
for (k=0;k<8;k++) { 


x1 = gen_xi( k, g->xi1 ); 
yl = gen _yi( k, g->yi ); 
1f( ( node[xi] [yi]->state != OBSTACLE) && 


( node[xi] [yi]->dist_from_start < best ) ) { 
best = node[xi] [yi] ->dist_from_start; 
best_ptr = node[xi] [yi]; 
} 


} 
g = best_ptr; 
g->state = X; 


116 


tdisplay.c 


/* 
File: taisplay.¢c 
Programmer: g.b. parker 
Environment: any 


Language: Cc 

Date: oe yuly $2 

Revised: 

Comments: Functions called by all searches to display the search on the 


IRIS. This file should not be excluded from Makefile if compiled on the SUN. 
ay 


#include “ga_search.h" 
#include <gl.h> 
#include <device.h> 


/* kaka K KK KKK KKK KK KK initialize KKK KKK KKK KKK KKK KKK KKK KKK KKK KK KKK KKK * / 


/* Initializes graphics systems for output */ 


initialize(title) 
eharee it le (33); 
{ 

/* set up a preferred size and location for the window */ 
prefsize (XMAXSCREEN+1, YMAXSCREEN+1-256) ; 
BPECEDOSItELON(O, 900, 0, 980) ; 

/* open a window for the program */ 
winopen ("search") ; 
/* put a title on the window */ 
wintitle(title) ; 
/* put the machine into double buffer mode */ 
doublebuffer () ; 
/* set RGB mode for color */ 
RGBmode () ; 
/* configure the IRIS (means use the above command settings) */ 
Geontig({); 
/* queue the redraw device */ 
qdevice (REDRAW) ; 
/* queue buttons needed */ 
qdevice(BUT6); /* ESC */ 
qdevice(BUTS50); /* enter */ 
qdevice(BUT4); /* right shift */ 
qdevice (BUT73); /* down arrow */ 


/* set the world coordinate system */ 
@rtho2 (-150766.0-—-1.0,66.0); 
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[* kkk KKK KK KKK KKK KKK draw terrain KR RRR RRR KKK KKK KKK KKK KKK KKK KK KK KK KK Ok / 


/* Called by searches to draw the node array */ 


draw _terrain(node, start,goal,current, dist, chrom) 
struct node rec *node loGiiioc), 
struct node rec ™™’ start, “goal, *Cunrent, 
float dist; /* dist_traveled */ 
struct factor struct chrom, 7 ind chrometaccome:, 
{ 

short value; 

Static int cont = T; 

int mmouse = F; 

inte first = 91, 

intwedo print; 


if( adjacent (start,current) || adjacent (goal,current) ) 
COnt) = 1; 
while( (mmouse || first) && cont ) { 


do print = F; 
draw grid(); 
draw_nodes (node, start, goal, current) ; 


while( qtest() ) 
Switch( qread(&value) ) 

{ 

case BUT6: /* “ECS" to terminate display for that search */ 
CONnts——5e, 
break; 

case BUTSO0: /* “return” to halt display us, 
mmouse = T; 
break; 

case BUT4: /* "shift" to continue display */ 
mmouse = F; 
break; 

case BUT73: /* “down arrow" to print node info to standard outpur 

tt f 

dou prinkg = 1; /* node is selected by mouse position */ 
break; 

default: 
break; 

} 

show_mouse (node, dist, chrom,do print); 


swapbuffers(); /* change the buffers ... */ 
first = F; 
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/* Kaa K KK KKK KKK KK show mouse KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK 


/* Shows mouse position and prints node info if selected */ 


show_mouse (node, dist,c,do_print) 

struct node_rec *node[66] [66]; 

eloat dist; 

Gmever factor Struct c; /* ind chrom_factor */ 
fat do print; 


{ 


int mx_ pix = getvaluator (MOUSEX) ; 
int my pix = getvaluator (MOUSEY) ; 


meemx = ((67 * mx pix)/980) - 1; 
imeemy — ((67 * my pix)/980) - 1 


° 
, 


if ( mx > 65 ) 


mx = 65; 
teemy > 65 ) 
my = 65; 


Reecolor (0,0, 0) ; 
eauare( node[mx] [my]->x, node({mx] [my]->y, 0.35 ); 


geedo. print ) { 

eernte("\n dist=%f", dist); 

printf ("\n(%d, %d) \ntd state\ntd btstate\ntf subtotal\ntf %d start 
\nt£f td goal\n%tf $d current\n%d frontier", 
node [mx] [my] ->xi, node[mx] [my]->yi, node[mx] [my] ->state, 
node [mx] [my] ->back_track state, node[mx] [my]->subtotal, 
node [mx] [my]--dist from start, c.start dist, 
node [mx] [my] ->dist_ from_goal, c.goal dist, 
node [mx] [my] ->dist_from_current, c.current_dist, 
node (mx] [my] ->frontier index); 

printf ("\ntf from below", node[mx] [my]->subtotal - 


aif 


(node [mx]/{my]|—->dist from start™ c.start dist + 


node [mx] [my] ->dist_from_goal * c.goal dist) 
printf ("\ncs=%d, cd=%d, ma=%d, m=%d\n", 
c.crowd sides,c.crowd diag,c.move_away,c.momentum) ; 
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es 


/* KKK KKK KEKE KKKKK draw nodes KK KI KKK KK KKK KK KKK KKK KKK KKK 


/* Draws node info, can be changed for color */ 


draw_nodes (node, start, goal, current) 
struct node_rec *node[66] [66]; 

struct nede rec “start, *Goal, “Current; 
{ 


aerite elegy les 


for (xi=0;x1i<=65; xitt) 
for(yi=0O;yi<=65;yitt) { 
Switch( node[xi] [yi]->state ) 
{ 
case OBSTACLE: 
RGBco lor (0, 0,0) ; 
squaref( node[xi] [yi]->x, node[xi] [yi]->y, 0.5 ); 
break; 
case VISITED: 
/* RGBcolor (0,0, 255))2 47 
circf( node[xi] [yi] ->x, node[xi] [yi]->y, 0.1 ); 
etre ( node[xi] [yil—2x, node l2uiiy i) =>77 Ce 500. 
break; 
case FRONTIER: 
/* RGBcolor (0,255,0); */ | 
circ( node[x1i] [yi]—>x, node[xi] [yil<=>y,70-3 ); 
break; 


/* RGBcollor (255, 07255). 7 
circf( current->x, current->y, 0.3 ); 


/* RGBeollor(2Z55, 0, 0s, 


circf( start->x, start->y, 0.4 ); 
circf( goal->x, goal->y, 0.4 ); 
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et! 





[III RII draw Grid 77 IOI III KK te / 


/* Draws the cross lines for the grid */ 


draw grid() 
{ 
int 1; 
elloat fi; 


/* draw the background color */ 
Pespeelor (255,255,255) ; 
clear(); 


RGBcolor(0,0,0); 


for (i=0;i<=66;i++) { 
fi = (float) (1-0.5); 
move2 (-0.5,f1); 
diraw2 (65.5, £1); 
moved (£1.,-0 .5) ; 
Grawz (£47 65.5)>; 


/* KRaAKKKKKKKKK KK KKK K squaref KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KK x / 
/* display filled square for 2D displays */ 

wold Ssquaref (xc, yc, d) 

moat xe, ye; /* center point of square */ 

Bloat d; /* half of side length */ 


{ 
rect f (xc-d, yc-d, xc+d, yc+d) ; 


/* KaKKKK KKK KKK KK KKK square KKK K KKK KKK KKK KKK KKK KK KKK KKK KKKK KKK KKK KKK x / 
/* display square for 2D displays */ 

void square (xc, yc,d) 

float xc,yc; /* center point of square */ 

Zloat .d; /* half of sidelength */ 


{ 
BeCt (XC-G, Vve-a, XCtG, YCtd) ; 
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tprint.c 


/* 
File: Epil ac 
Programmer: g.b. parker 
Environment: any 


Language: e 

Date: 9 july 92 

Revised: 

Comments: Prints to standard output 


= 
#include "ga_search.h" 


/* REKKKKKKKKKKKRKRKRKK print density KKK KKKKKEKKKEKKKEKEKKEKKKEKEKKEKKKKKKKKKKK x / 


/* Prints density terrain to standard output */ 


print density( density ) 
int density([(16] [16]; 


{ 
int? i; sar 


for ( j=15; j>=0; j-- ) {° 
printf (7 \ns, 
for ( 1=0734<=157 911 te) 
Printe ("Sx “density | iiigi)- 


/* KKK KKKKK KKK KK KKK print population KRHaKEKEKKKKKKKKAKKKKKKKKKKKK KKK KKEK x / 


/* Prints the population to standard output */ 


print _population( i ) 
SCrIUCE iIndiviaual Strucem. scl, 
{ 

int 


for (k=0;7k<32; kt+) 
printf ("\n td x td t£ td", k, if({k]->chrom.alleles, i[(k]->fitness, i[k]- 
>fit_sum, i[(k]->previous_index) ; 
} 
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/* KRaKKK KK KK KKK KR KKK Drine node kkk kak Ke KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK K « / 


/* Prints node terrain to standard output */ 


print_node( node ) 
Struct node_rec *node[66] [66]; 
{ 


mit iL, J; 


for ( j=65; j>=0; j-- ) { 
eronte (™\n") ; 
moo ( 1=0; 1<=65; it+ ) 
Switch (node[i] [j]->state) 

{ 

case UNTOUCHED: 
PeIneet.,” ):; 
break; 

case OBSTACLE: 
prinerc' #" ) > 
break; 

case VISITED: 
Bein e6" oO"); 
break; 

case FRONTIER: 
Primer ( f).; 
break; 

case START: 
Brie (os. ) 5 
break; 

case GOAL: 
Prince (GCG); 
break; 

case CURRENT: 
DPeinte (©) ; 
break; 

case SHORTEST: 
DYINEE Gs) ; 
break; 

case X: 
Demme r(x); 
break; 
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